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I .   INTRODUCTION 

Satellite  communications,  especially  digital,  are  well 
under  way.   One  of  the  many  observations  that  satellites 
provide  are  the  meteorological.   These  give  the  location  of 
a  storm  or  a  typhoon  in  terms  of  geological  coordinates.   It 
is  possible  to  have  images  of  the  store  or  typhoon  cloud 
cover . 

These  observations  occur  quite  randomly  in  time. 
Meteorologists  like  to  present  their  forecasts  equally 
spaced  in  time. 

The  attempt,  here,  is  to  try  to  compare  an  adaptive 
Kalman  filtering  algorithm  to  estimate  the  location  of  the 
eye  of  a  storm  with  optimum  track  values  that  are  given  from 
the  meteorologist's  analysis.   The  filtering  process  combines 
all  the  available  measurement  data  with  prior  knowledge  about 
the  system.   It  produces  an  estimate  of  the  location  of  the 
storm  (latitude,  longitude)  in  such  a  manner  that  the  mean 
square  error  is  minimized.   The  parameters,  which  are  essential 
design  elements  of  a  Kalman  filter,  are  the  measurement 
noise  covariance  matrix,  R,  the  excitation  covariance,  Q, 
the  initialization  covariance  of  error  in  the  filter  itself, 
P(l/0),  and  the  transition  matrix,  0. 

In  most  physical  processes  that  one  desired  to  track, 
many  of  these  parameters  change  during  the  process  of 


tracking.   The  measurement  noise  associated  with  the 
observations  can  change  if  a  different  sensor  is  used,  or 
a  similar  sensor  obtaining  measurements  from  some  different 
geometry  relative  to  the  object  being  tracked.   If  the  object 
being  observed  is  acted  upon  by  external  forces,  then  the  Q 
matrix  should  be  changed  to  account  for  these  external 
excitations.   Most  processes  being  tracked  will  change  in  their 
dynamic  characteristics  during  the  observation  time  and  so  the 
transition  matrix  ideally  should  also  be  changed.   Further, 
the  time  between  observations  quite  often  occurs  randomly 
in  time.   All  of  these  things  bring  about  the  need  to  change 
the  Kalman  filter  to  adapt  as  the  process  changes. 

It  is  possible  to  change  the  parameters  of  the  Kalman 
filter  by  sensing  the  error  between  the  observation  and  the 
prediction  from  the  track  of  what  that  observation  should  be . 
If  a  gate  is  established  representing  95%  of  the  normal 
random  perturbations  of  the  process,  then  when  this  error 
exceeds  the  magnitude  of  the  gate,  one  can  reasonably 
ascertain  that  the  filter  is  no  longer  properly  representing 
the  observed  process.   In  the  work  attacked  here,  real  data 
was  obtained  from  satellite  observations  and  the  qualitative 
observation  error  was  established  (PCN  #). 

The  error  covariance  matrix  in  terms  of  error  ellipsoids 
along  the  track  gives  a  measurement  for  the  worthiness  of 
the  algorithm  parameters. 


II.   KALMAN  FILTERING  TECHNIQUES 

In  a  linear,  discrete  system,  the  state  and  measurement 
equations  are  given  by 

x(k+l)  =  0(k)x(k)  +  £w(k) 

and 

z(k)  =  H(k)x(k)  +  v(k) 

where  x  is  the  state;  $   is  the  transition  matrix;  _T  is 
excitation  noise  matrix;  H  is  the  measurement  matrix;  w 
and  v  are  the  excitation  and  measurement  noise  correspondingly, 
assumed  uncorrelated,  zero  mean  white  Gaussian: 

E[w(k)  -wT(j)]  =  Q(k)6kj 
and 

E[v(k)  -vT(j)]  =  R(k)5kj 
and 

E[w(k)  ]  =  0,  E[v(k)  ]  =  0 


where  Q(k)  and  R(k)  are  covariances  of  excitation  and 
measurement  noise.   Now  if  x(k)  is  the  estimated  state 
value  after  the  kth  measurement  and  x(k|k-l)  is  the 
predicted  value  of  the  state  before  the  kth  measurement 
we  have : 


x(kjk-l)  =  f (x (k-l| k-1)  ,k-l)  ,  where  f  is  any  function 


The  state  error  vector  is  defined  to  be 


x  (k)  =  x  (k)  -  x  (k) 


and  the  predicted  state  error  vector  is  defined  to  be 


x*(k|k-l)  =  x(klk-l)  -  x(k) 


The  covariance  of  state  error  matrix  is  defined  to  be 

P(k|k)  =  E[x(k)  -xT(k)  ]  , 
and  the  predicted  covariance  of  state  error  is  defined  as 

P(k|k-1)  =  E['x(k|k-1)  -xT(k|k-l)  ]  . 
The  state  excitation  matrix  is  defined  by 

Q(k)  =  F(k)E[w(k)  -wT(k)]  -rT(k)  , 
and  the  measurement  noise  covariance  matrix  is  defined  by 

R(k)  =  E[v(k)  -vT(k)  ]  . 
The  Kalman  filter  equations  are: 

P(k  +  l|k)  =  0(k)P(k|k)0T(k)  +Q(k) 

G(k)  =  P(k|k-l)HT(k) [H(k)P(k|k-l)HT(k)+R(k) ]_1 

P(klk)  =  [I-G(k)H(k)]P  (k|k-l) 


10 


x(k|  k-1)  =  0(k)x(k-l|  k-1) 

| (k|k-l)  =  H(k)x(k|  k-1) 

x(k|k)  =x(k|k-l)  +  G(k)  [z  (k)-z  (k|k-l)  ] 

The  initial  condition  of  P  (error  covariance  matrix)  and  the 
Q  and  R  matrices  are  the  determining  factors  in  the  filter 
structure.   For  Q  having  main  diagonal  values  greater  than 
R  means  that  we  have  greater  uncertainty  in  the  state 
estimate  than  in  the  observation.   Thus  the  new  state 
estimate  is  more  dependent  upon  the  new  measurement  and  less 
related  to  prior  estimates.   The  inverse  is  also  true.   For 
R  having  greater  diagonal  terms  indicates  that  the  new 
measurement  are  subjected  to  stronger  corruptive  noises, 
and  so  should  be  weighted  less  by  the  filter.   The  gains  (G) 
are  lower.   The  P  is  responsible  for  the  initial  transient 
performance  of  the  filter. 
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III.   ERROR  ELLIPSOIDS 

The  error  covariance  matrix  in  each  stage  of  a  Kalman 
filter  process  gives  insight  into  the  quality  of  the  track 
occuring . 

The  diagonal  terms  (P,  ,  and  P-,-?)  are  the  variances  of 
uncertainty  in  our  knowledge  of  latitude  and  longitude.   Their 
respective  off  diagonal  terms  are  covariance  between 
latitude  and  longitude. 

The  square  roots  of  the  diagonal  terms  gives  us  the 
rms  errors  in  out  estimates  of  longitude  and  latitude.   Having 
the  definition  of  the  structure  we  are  dealing  with  (in  our 
case  the  satellite  observations)  and  its  uncertainties 
(expressed  by  the  PCN  number-actually  by  the  values  of  R) 
we  can  see  how  the  K.F  performs  through  its  error  covariance 
matrix.   Expressing  the  P  matrix  in  an  ellipsoid  of  constant 
probability,  one  obtains  a  visual  appreciation  for  the 
worthiness  of  the  algorithm  parameters.   The  representation 
requires  that  the  errors  are  normally  distributed. 

The  joint  probability  density  function  is: 

e-l/2eT(ki  k-l)P_1  (k|k-l)e(k| k-1) 

where  e(k|k-l)  is  the  predicted  state  error  vector.   Setting 
the  exponent  equal  to  a  constant  value,  we  are  going 
to  have  a  curve  which  is  an  ellipse.   This  ellipse, 
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however,  does  not  have  its  major  and  minor  axis  aligned  with 
the  coordinate  system.   Instead  its  axis  (x1,  y1)/  comes 
from  the  following  transformation: 

x'  =  xcos0  +  ysintf 

y'  =  ycostf  -  xsin# 

where  x,  y  are  the  latitude  and  longitude  in  our  case  with 

6  =    1/2    tan"1!    2^OV(x-y)] 

(7        -     <T 

x  y 


where   cos(x,y)    =   P    „ (k   k-1) 


C72    =   P..  (klk-1)    and   cr2    =    P.-(klk-l) 
x  11       I  y  22       1 


The  new  variances  are  calculated  by 


2     2 
,  2     x     y    cov (xy ) 
a  x  2  sin2  0 


rr2  2 

, 2     x     y    cov (xy ) 
°"y      2       sin22 


Incorporating  the  above  equations,  error  ellipses  are 
presented  in  subsequent  figures  with  the  satellite  tracks 
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IV.   SATELLITE  TRACKING-SIMULATION  RESULTS 

Data  from  the  Annual  Tropical  Cyclone  Data  for  the 
Typhoon-Nelson  appear  in  Tables  1  and  2  in  terms  of 
longitude  and  latitude.   The  best  track  data  appears  in 
six-hour  intervals  for  twelve  days  and  the  satellite  fixes 
(observed)  in  random  time  intervals. 

The  data  in  the  K.F  algorithm  parameters  are: 
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X(l/0)  = 


In  Figure  1  we  have  a  representation  of  the  best  track 
observed  and  the  K.F  track.   An  error  ellipsoid  at  the  25th 
stage  of  the  process  appears.   It  seems  that  K.F  follows 
the  observed  track  more  closely  than  the  best  track.   This 
is  due  to  the  adaptive  gating  and  Q  relative  to  R. 

The  gain  history  for  0(1,1)  fluctuates  between  0.7-0.8, 
never  arriving  at  a  stable  value.   The  G(3,l)  reaches  a 
stable  value  of  0.12. 

The  above  appears  in  Figure  2 .   The  track  for  prediction 
appears  in  Figure  3 . 

The  errors  EB1  and  E0B1  represent  (YH-BLAT)  and 
(LAT-BLAT)  respectively.   They,  along  with  EB2  and  E0B2 
appear  in  Table  3  and  Figures  4  and  5. 

Implementing  Julian  Time,  we  have  a  comparison  in  common 
time  for  24  points  only.   It  can  be  seen  that  in  terms 
of  the  latitude  error  the  K.F  is  close  to  the  best  track 
values . 

In  an  attempt  for  better  performance  the  Q 
matrix  was  changed  in  the  algorithm.   Also  the  innovation 
errors,  in  terms  of  latitude  and  longitude,  exceeding  the 
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magnitude  of  the  gate  (v/P(k  k-1)  +  R)  resulted  in  changing 
the  values  of  the  gains.   The  above  correction  makes  the 
filter  more  adaptive  now  and  the  errors  EB1  and  EB2  appear 
smaller  on  an  average.   The  representation  of  the  above 
correction  in  terms  of  trajectories,  gains,  predictions  and 
errors  appear  in  Figures  6,  7,  8,  9,  and  10  and  Table  4. 
The  computer  program  appears  in  Appendix  A. 
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TABLE  1 
OBSERVED  DATA  (SATELLITE  FIXES) 
TIME  LAT         LONG 


1804. 00 

3.70 

160. 90 

1809.00 

4.00 

160.  10 

1318.00 

4.70 

157.  70 

1900.00 

5.30 

155.50 

1903.48 

5.  80 

154.30 

1906. 00 

5.80 

153. 90 

1909- 00 

6.  00 

152.60 

1916.33 

6.  70 

150. 70 

1921 .00 

6.80 

150. 00 

2000. 00 

7.20 

149. 60 

2003. 00 

7.20 

149.30 

2005. 18 

7.50 

149. 10 

201 2. 00 

7.00 

147. 40 

2016.21 

7.20 

146. 10 

2018. 00 

7.20 

146- 00 

2021 .00 

7.20 

145. 70 

2100. 00 

7.30 

145.20 

2103. 00 

7.50 

145. 00 

2105- 06 

3.20 

1 44.40 

2106. 00 

8.50 

144. 00 

211 2.00 

7.80 

143. 00 

21 16. 00 

3.  00 

14 1. 90 

211 7.51 

8.  00 

141. 70 

2200. 00 

8.30 

139.30 

2203.00 

8.  50 

1 39.00 

220  4.54 

8.50 

138. 10 

2206.00 

fa.  t>0 

138.00 

2212. 00 

8.60 

1 35.90 

2^1 6.  00 

8.30 

134.70 

221 7.40 

9.00 

134.20 

2300. 00 

8.80 

1 33.60 

2306. 00 

9.30 

132.30 

2312. 00 

9.40 

131.40 

231 7.23 

9.50 

130.  10 

2321 .00 

9.50 

129.60 

2400.00 

9.50 

129.50 

2403.00 

9.oO 

129.20 

2406. 13 

9.80 

128.70 

2412.00 

9.70 

128.20 

241 6- 00 

9.60 

127.30 

2413.00 

9.  70 

127.60 

2421.00 

9.90 

127.40 

2500. 00 

10.00 

127.30 

2503.00 

10.  10 

127. 10 

2506.01 

10.  1  0 

126-30 

2509. 00 

10.30 

126.70 

2512.00 

10.30 

126.30 

2518. 00 

10.40 

125.50 

2521 .00 

10.40 

125.  10 

2600.00 

10.30 

124.50 
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TABLE  2 

BEST 

TRACK  DATA 

TIME 

3LAT 

BLONG 

I  806*  00 

3.80 

160. 70 

1312.00 

4.50 

158.80 

181 8. 00 

4.90 

157.  10 

1900.00 

5.50 

155.50 

1906. 00 

5.90 

153.90 

1912.00 

6.40 

152.40 

1918. 00 

6.70 

151-00 

2000. 00 

6.90 

149.90 

2006. 00 

6.90 

143. 70 

2012.00 

7.00 

14  7.  60 

2018. 00 

7.30 

14.6-50 

2100. 00 

7.50 

1*5.50 

2106. 00 

7-80 

144.40 

2112.00 

7.90 

142.90 

2118. 00 

a.  10 

141 .30 

2200.00 

8.30 

139. 50 

2206.00 

8.60 

137.70 

2212.00 

8.80 

136. 10 

2213.00 

9.00 

134..  ao 

2300-00 

9.10 

133.50 

2306. 00 

9-20 

132-40 

2312.00 

9.30 

131.40 

231 6. 00 

9.50 

130-40 

2*00.00 

9.  70 

129.50 

2*06. 00 

9.30 

128.80 

2412.00 

9.80 

128-20 

24.13.00 

9.90 

127.70 

2500. 00 

10-10 

127.20 

2506. 00 

10.30 

126.30 

2512.00 

10.50 

126-20 

251 8. 00 

10.40 

125.50 

2600. 00 

10.30 

124.80 

2606. 00 

10.50 

124.00 

2612. 00 

11  .00 

123-50 

2618. 00 

11  .50 

122.90 

2700. 00 

1  1  .70 

122.20 

2706. 00 

11.30 

121-30 

2712. 00 

1  1  .90 

120-40 

2718.00 

12.00 

119.50 

2800. 00 

12.40 

11 8.50 

2606- 00 

12.90 

11  7-80 

2812.00 

13-30 

11  7-20 

2813. 00 

13.70 

11 6-oO 

2900. 00 

14.  10 

116. 10 

2906. 00 

14,30 

11 5.80 

2912. 00 

14.-20 

115.40 

2918. 00 

14.00 

1 1 4- 90 

3000. 00 

13.90 

11 4.50 

3006. 00 

13.30 

1 1 4.20 

3012- 00 

13.50 

1 14. 00 
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TABLE  3 
ERROR  VALUES  IN  TERMS  OF  LATITUDE  AND  LONGITUDE 
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TABLE  4 
ERROR  VALUES  IN  TERMS  OF  LATITUDE  AND  LONGITUDE 


JTT    M" 

-  -  1 

if  32 

rr 

1 

erJ?2 

7  20  j4  2 

-0.  Jl 

0  .  ol 

— l 

.20 

0  .  60 

72'in±<l 

-j  .27 

0  .  0 .5 

-  J 

.^0 

0  .  2 

~"  ?  tj  o '"  4 

0.16 

-0.12 

-0 

.  10 

0  .  0 

72p o72 

0.26 

-0.57 

r\ 

.  30 

-0  .  30 

720&3* 

0.  13 

-0.21 

j 

.  Li 

-0.20 

7^0  ^  30 

-0.12 

—  '     .  o  a 

—  j 

.  i  j 

-0  .50 

7^.06^6 

-0  ,  .2  1 

-  7 . T  ft 

— ) 

.  20 

-0.30 

720  70  2 

0  .  7  * 

-0  ,~7 

:) 

.  7  1 

-0  .  4,] 

7L071B 

0.  10 

■">  .J  J 

—  j 

.  1° 

~>  .  1  0 

7^7^ 

-0.0  '" 

-  0  .  1  c 

; 

•  vj 

-0.2  0 

7  20  7."!  j 

-  o  .  1 6 

0  .  ~j 

—  j 

.   10 

0  .  3u 

72°  71  ^ 

-  0  .  o  0 

-0.13 

—  ■) 

.  2  0 

-0  .  2  J 

7  c  }  7a4 

-0.17 

0  .  I  V 

— 0 

.  30 

0  .  10 

7^0  ^0 

0  .  1  o 

-  0  .  0  ft 

w 

.  1  1 

-o.  to 

7_  )  7  "n 

0.  1  d 

-0.04 

3 

.  1  "> 

0  .  0 

72f>  70d 

-  0  .  14 

-0.01 

—  J 

.  20 

0  .  0 

7  2  0  7  7  A 

- 0 .  on 

-  o  .  i  o 

) 

.  0 

-0.10 

7.;  17*0 

-0.  04 

o.o^ 

—  -j 

.  1  0 

r>         rt 

7  2  0  7  3a 

—  u  .  ^G 

-0.1  I 

— ' 

-olio 

720 70  2 

-0.21 

-  0  .  0  7 

— fj 

.  in 

0.  1  j 

'20  7")  -j 

-0.21 

-  0  .  C  9 

-  .. 

.20 

o .  0 

720  30«t 

-0.  io 

-0.06 

—  ) 

.  2n 

0.  lu 

7cOcllu 

-0.02 

0.09 

i 

.  0 

0  .0 

7^0^! ^ 

-0.00 

-"  .10 

i 

.  0 

-0.30 

28 


LJ 


E-h 


w 

O 

w 

P 
+J 

■H 
-P 
rd 


0) 
S-l 

3 

•H 


1903^93 


29 


00 

o 

CM 
IN, 


CD 
CO 
[N 

o 
cm 

N 


CM 
CO 

in 


i r 

Z        9"I       Z'\ 


1 1 

CM 

m 

IX 

5-1 

0 

u 

u 

CO 

H 

to 

IS. 

QJ 

o 

T3 

™   r   i 
in    Lu 

z: 

3 
+J 
•H 

1 — ' 

Cn 

C 
0 

CN 

J 

o 

CM 

Cn 

o 

I-l 

a 

o 

>-i 

en 

CO 

o 

3 

CM 

•H 

x 

fa 

CO 

CD 

CO 

a 

CM 

IN. 

CM 

t1 

CO 

o 

CM 

IN 

VQ  0        fr'O-     8"0- 

sao3'zg3 


30 


V.   RANDOM  TRACKING 

To  examine  the  adaptability  of  the  K.F,  another  storm  was 
created.   The  equations  for  this  simulated  storm  were: 


BLAT(s)  =  BLAT(s-l)-V  -T 

2C 


BLONG(S)  =  BL0NG(S-1)-V  -T 

y 


LAT ( S )  =  BLAT ( s )  +  V ( s ) 
LONG(S)  =  BLONG(s)  +  V(s) 

where  T=6hr,  V  =10°/24hr,  V  =5°/24hr  and  V  =  measurement 

x  y 

noise  (created  by  a  random  generator  subroutine) . 

Implementing  the  above  equations,  a  "true"  and  an 
"observed"  track  were  created.   These  two  data  files  were 
named  BESTRACK  and  OBSERVED  and  appear  in  Table  5  and  6 . 

Simulating  with  the  above  new  data,  the  K.F  algorithm 
appeared  to  track  well.   Figure  11  shows  the  "true", 
"observed"  and  K.F  tracks. 

Figure  12  indicates  the  gain  history  in  terms  of  Gil 
and  G13.   This  approach  gave  stable  values  of  0.27  and 
0.01  respectively  after  the  4th  discrete  point  in  time. 
After  this  time  the  gain  does  not  vary  any  more.   This 

means  that  the  innovation  error,  (z (k) -x (k| k-1)  )  is 

i 
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weighted  each  time  by  the  same  quantity  after  the  4th 
observation.   Having  lower  values  in  the  diagronal  terms  of 
the  Q  matrix  in  comparison  with  R,  in  this  case,  means  that 
we  have  greater  uncertainty  in  the  measurements  (observed 
data)  relative  to  the  model  uncertainties.   So  the  gains  are 
smaller  and  the  filter  no  longer  "tracks"  the  measurements 
closely . 

As  far  as  the  latitude  and  longitude  errors  are  concerned 
it  can  be  seen  that  EB1  (YH-BLAT)  and  EB2 (XH-BLONG)  are 
nearly  zero  and  are  smaller  in  comparison  with  E0B1  (LAT-BLAT) , 
E0B2  (LONG-BLONG) .   This  shows  the  ability  of  the  algorithm 
to  follow  the  "true"  values  more  than  the  measurement  if  the 
latter  has  been  corrupted  with  noise. 

The  error  values  and  the  plots  appear  in  Table  7  and 
Figures  13  and  14. 
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TABLE  5 
BESTRACK  DATA 

TIME  SLAT  BLONG 

1900.00  160,00  3,30 

1806-00  157.50  2.80 

1812.00  155.00  1.80 

1813.00  152.50  0.80 

1900.00  150.00  -0.20 

1906.00  147.50  -1.20 

1912.00  145.00  -2.20 

1913.00  142.50  -3.20 

2000.00  140.00  -4.20 

2006-00  137.50  -5.20 

2012.00  135.00  -6.20 

2013.00  132.50  -7.20 

2100.00  130.00  -8.20 

2106.00  127.50  -9.20 

2112.00  125.00  -10.20 

2119-00  122.50  -11.20 

2200.00  120.00  -12.20 

2206.00  117.50  -13.20 

2212.00  115.00  -14.20 

2213.00  112.5  0  -15.20 

2300-00  110.00  -16.20 

2306.00  107.50  -17.20 

2312.00  105-00  -18.20 

2318.00  102.50  -19.20 

2400.00  100.00  -20.20 

2*06-00  97.50  -21.20 

2412.00  95.00  -22.20 

2413.00  92.50  -23.20 

2500.00  90.00  -24.20 

2506.00  37.50  -25.20 

2312.00  85.00  -26.20 

2518.00  82.50  -27.20 

2600.00  30.00  -29.20 

2606.00  77.50  -29.20 

2512.00  75.00  -30.20 

2613.00  72.50  -31.20 

2700-00  70.00  -32.20 

2706.00  67.50  -33.20 

2712.00  65.00  -34.20 

2713.00  62.50  -35.20 

2800.00  60.00  -36.20 

2305-00  57.50  -37.20 

2312.00  55.00  -38.20 

2818.00  52.50  -39.20 

2900.00  50.00  -40-20 

2906-00  47.50  -41.20 

2912.00  45.00  -42.20 

2913.00  42.50  -43.20 

3000.00  40.00  -44.20 

3006.00                  37.50  -45.20 
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TABLE  6 
OBSERVED  DATA 

TI*E  LAT  LONG 


1800,00  160.00  3.30 

1306,00  157.30  2.60 

1812.00  15*. 92  1.72 

1815.00  152.71  1.01 

1900.00  150.03  -0.17 

1906.00  14.7.63  -1.07 

1912.00  14.4. 95  -2.25 

1918.00  142.61  -3.09 

2000.00  139.99  -4..  21 

2006.00  137.26  -5.44 

2012.00  134. 6.2  -6.58 

2018.00  132.39  -7.31 

2100.00  130.31  -7.89 

2106.00  127.37  -9.33 

2112.00  125-19  -10.01 

2113.00  122.44  -11.26 

2200.00  120.19  -12.01 

2206.00  117.40  -13.30 

2212.00  114.87  -14.33 

2218.00  112.92  -14.78 

2300.00  110.03  -16.17 

2306-00  107.65  -17.05 

2312.00  104.83  -18.37 

2318.00  102.96  -18.74 

2*00.00  99.90  -20.30 

2406.00  97.65  -21.05 

2*12.00  95.29  -21.91 

2418.00  92.44  -23.26 

2500.00  90.51  -23.69 

2506-00  87.13  -25.52 

251*;. 00  84.84  -26.36 

2518.00  82.40  -27.30 

2600.00  30.00  -28.20 

2606.00  77.78  -28. 9<; 

2512.00  75.08  -30.12 

2613.00  72.50  -31.20 

2700.00  69.96  -32.2a. 

2706.00  67.22  -33.48 

2712.00  64.79  -34.41 

2713.00  62.75  -34.95 

2800.00  59.33  -36.37 

2305.00  57.37  -37.13 

2812.00  54.72  -33.48 

2318.00  52.43  -39.27 

2*300.00  50.06  -40.14 

2905.00  47-91  -40.79 

2912.00  44.33  -42.32 

2913.00  42.53  -43.12 

3000.00  40.14  -44.06 

3006-00  37.2b  -45.44 
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TABLE  7 
ERROR  VALUES  IN  TERMS  OF  LATITUDE  AND  LONGITUDE 


J  1  I  ;M= 

FBI 

c*2 

eoj  l 

F  U  8  2 

7  20624 

-0.  OS 

-0  .00 

J  .0 

0  .0 

7  20  6  3  0 
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0.31 

7^0  70  2 

-  0  .  05 

-O.Oo 

-U.M 

-0.13 

720 70 a 

0.  02 

o.oi 

0.1^ 

0  .  1  ; 

7^371 4 

0.00 

-o.oo 

-0  .05 

-0.05 

720  ^0 

0.  0  6 

O.Oj 

0.10 

0  .  19 

720 72 o 

0  .  02 

0  .0^ 

-0.  1  0 

-0.10 

7  20  7^2 

-0.02 

-0.02 

-0.13 

-0.13 

^207Jrt 

0.  10 

0.  10 

J  .  4  2 

0.42 

720 744 

O  .  10 

0.0S 

0.03 

0.03 

720  75(J 

0.12 

0.12 

u  •  1  ~ 

0.16 

7  20  75  6 

0.05 

0.05 

-3.17 

-0.17 

7  2  0  7  f .  2 

0.17 

0.17 

J  .  4  6 

0  .46 

720 7Sd 

0.  1  1 

0  .  1  l 

-0  .  1  0 

,      -  0  .  10 

7^0  774 

0.   12 

0.  1  <L 

0.  15 

0.  1  j 

7  20  78  0 

O.  1  & 

0.1S 

j  .29 

0  .  ?9 

720  78o 

0.12 

0.12 

-i  J  .05 

-O.OO 

7 JO  70  2 

0-23 

0  .23 

0.51 

0  .51 

7^07Qd 

0.  10 

0.10 

-J  .  32 

-0  .32 

7^0804 

0.  02 

0.0  2 

-0.16 

-0.16 

72  03  10 

-0.02 

-0.02 

-0  .  in 

-  0  .  10 

7^061 o 

-o.o^ 

-0.03 

0.0 

0  .0 

r^0^22 

0  .  04 

0.04 

■j  «  2  3 

0.28 

72  0  82  ti 

0.  0  5 

0.05 

3  .  L,  J 

0  .Oo 

720  -j  3^ 

0.  04 

0.0  4 

0  .  0 

0  .0 

7^  0  2>40 

0.01 

0.01 

-0  .04 

-0.04 

7^3-i4  6 

-0.07 

-0.07 

-0  .28 

-o  .  2a 

720  do 2 

-0.  12 

-0.12 

-J.il 

-0.21 

7  20 658 

-0.04 

-0.04 

j  .25 

^.2  3 

7^:0  ri(i<t 

-0.  08 

-O.Ob 

-O.I7 

-0.17 

7 2  Oc.  7 0 

-0.05 

-0.05 

0  .07 

0.0/ 

72  0  676 

-0.  1  2 

-0.12 

-  )  .28 

-0.26 

72  OciP.  2 

-0.  1  1 

-0.  1  1 

-U   .07 

-0  .0  7 

72  0*48  ti 

-0.08 

-  0  .  0  8 

'J  .  0  6 

0.0  5 

7^.0  rf  9  4 

0.  05 

0  .  0  b 

0-4  1 

0.41 

7^0900 

0.  0  2 

0.0  2 

-0.12 

-0.12 

7  2  0  9  0  o 

0.04 

0.04 

0  .  C  8 

0.08 

7^0912 

0.  OS 

0.0  3 

0.14 

0.14 

7  j  0  V  1  8 

-  0  .  IB 

-o.i  a 

-O   .24 

-0.2* 
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TABLE  7 
ERROR  VALUES  IN  TERMS  OF  LATITUDE  AND  LONGITUDE 


J  1  I  ,.ll 

ti»i 

c  ^  2 

cOO 

FU  M2 

7  J  n  o  ?  4 

-U.l)H 

-0  .  Oc 

J  .  o 

0  .  0 

7  2Mb  3  0 

-0.20 

-0.2  J 

-.)  .  <^o 

-0.20 

'^b  16 

-  0  .  1  2 

-  0  .  1  -J 

-  .  .  0  l1 

-0  .  OS 

7  J')o4^ 

0.  I  1 

0.10 

D  .21 

0.21 

7^0t,'id 

0.10 

0  .  o  e, 

o.oi 

0  .Oj 

7^  Oof.** 

o.  ir- 

0.  1  J 

0  .    1    1 

0.  1   « 

720^60 

0.  0') 

0.07 

-J  .Ori 

-O  .0  3 

7  jOf,of> 

0.   1  1 

0  .  !  0 

I) .  1  I 

0.  11 

7JI!67<: 

0.  0  J 

0.0  7 

-0.01 

-0.^  1 

7^0  u7e5 

-0.  02 

-  0  .  0  j 

-  )   .  2  4 

-  O  .  2  •♦ 

72  0  oH4 

-0.14 

-  o  .  i  --, 

-:  .  3r* 

-0  .  1   I 

rj'^o^J 

-  0  .   15 

-o.io 

-o  .  1  1 

-0.11 

7«i0vj'*b 

-0.  02 

-0.03 

0  .  J  t 

o .  u 

7  <;  0  70^ 

-0.0^ 

-O.Oo 

-.,11 

-0.1    J 

7  2  0  7  0  « 

o  .  o? 

n.o  i 

0  .  i  o 

0.  1   ) 

^in-. 

0  .  0  0 

-o.oo 

-0  .  0  ~> 

-o.oo 

7  20 7v 0 

0  .  lm  > 

n.o-j 

0  .   1  :> 

n.  w 

7  20 7 2o 

O  .  02 

o  .  o^ 

-v)   .    1    0 

-0.10 

7  2 0  71^ 

-o .  y 

-0.02 

-0.1"? 

-0.1    J 

72o  7i(i 

0.10 

0.10 

J  .  4  ?. 

0.42 

7^1  76  4 

o.io 

0.0* 

0.01 

0.0.1 

7  2  0  7SQ 

0.12 

0.12 

U   .   1    " 

0.  is 

720  7'30 

0.  OS 

0.05 

-  J .  1 T 

-0.17 

72  0 7? 2 

0.  1  7 

0.17 

>.<.'> 

n  .  4 1. 

720  7'«o 

0.11 

0  .  1  l 

-0.10 

-  0  .  10 

^077^ 

0.   1  2 

0.\     c 

0  •  l'"> 

O.   1    > 

7 2 17«0 

O  .  1  H 

o .  t  •; 

.  .  sn 

)  .  ?•) 

7  J   )  Mo 

0.12 

o.i2 

-•  J  .OS 

-  0  .  0  o 

72  0  7')^ 

U.2.7 

0  .    'J 

) .  r.j  l 

o  .  -  1 

7^0  7<»<; 

0.10 

0.10 

-  j  .  \1 

-0  ..12 

7^0  -in-'* 

O.  0? 

0.02 

-0  .  1  *> 

-T.h 

7  2  0  -j  I  0 

-0.02 

-0.02 

-J  .  1  n 

-  r»  .  1  0 

7_0Hlo 

-  0  .  0  7 

-  0  .  0  J 

i)  .  0 

0  .  0 

7<:0fV2 

O.ui 

v)  .04 

-<     •   <rj  " 

0  .  ?  M 

720  kJ2'l 

o.o*". 

0.05 

J  .  o  ■' 

0  .  0  o 

7  20    J  U 

0.0<* 

0.04 

Ci  .  0 

0  .  o 

7ijnj4i) 

0  .  O  1 

0.0  1 

- .  !   .04 

-O  .04 

"*  2  3  -\  '-•  b 

-0.07 

-  0  .  )  7 

-  "   .  2  * 

-  0  .   >  A 

7  2  1  ci  ">  2 

-  0  .   1  ? 

-o.  1  J 

-J  .  ^  I 

-0.21 

7  ,»  0  o  "i  M 

-O  .  04 

-0.04 

J   .  2<S 

r» .?  i 

7^  3  ri»S  <♦ 

-0.  OH 

-0.  Ob 

-j.r 

-0.17 

7^0o70 

-  o .  o  s 

_o  .  o  , 

0  .  0  7 

0.0/ 

7J0-j7o 

-  0  .   1  2 

-0.12 

-  )    .  2  * 

-0  .;»•} 

72  3<JS2 

-0.  1  1 

-0.  1   I 

-.1    .  O  7 

-  0  .07 

72Tl-?M 

-  0  .  0  s 

-  0  .  ( 1 1 ; 

.j  .06 

0.0  J 

T^n.fN 

o .  o  r> 

0  .Os 

0  .  -  1 

0.41 

7^0    10  0 

O.  02 

O  .  0  2 

-0.1"! 

-o.u 

7  i!0  O0o 

J.O'* 

O.0  4 

0  .OS 

o  .0  -3 

7^0   il  2 

0.0'? 

P  .  O.J 

0   .  1  4 

0.1* 

7  ^r  0  S>  l  M 

-  0  .  1  n 

-0.1  .0 

—  .'      .    c    & 

-0.2* 
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VI.   CONCLUSIONS 

The  K.F  approach  to  estimate  the  storm's  location  appears 
to  be  very  accurate.   That  comes  from  the  comparison  with 
the  meteorologist's  analysis  results.   Concerning  also  the 
fact  that  the  latter  was  performed  after  the  storm's  occurence 
one  can  see  the  advantages  of  the  K.F  algorithm.   During 
operation  of  the  filter  the  actual  residual  sequence 


[z_(k) -jz  (k|  k-1)  ]  is  compared  to  the  gate  [  v  P  (k  k-1) +R  (k)  ] 
which  actually  is  the  square  root  of  the  residual  covariance. 
Being  a  white  Gaussian  sequence,  the  residual  is  bounded  by 
this  gate. 

When  the  gate  is  exceeded  the  model  is  determined  invalid 
within  the  filter  and  a  modification  in  terms  of  Q  and  G 
takes  place  to  adapt  to  the  situation.   At  this  point,  and  if 
the  excess  occurs  in  only  one  component  of  the  vector  residual 
process,  one  can  further  deduce  that  the  measuring  device 
generating  the  particular  component  is  the  source  of 
difficulty  (a  sensor  failure)  . 

The  error  ellipsoids  of  the  process  also  give  insight 
into  the  filter  performance  in  a  more  general  case,  referring 
to  many  sensors  with  a  greater  variety  of  uncertainties,  the 
adaptive  K.F  algorithm  could  be  a  very  advantageous  approach. 
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APPENDIX  A 

COMPUTES  ALGORITHM 

KALMAN  FILIE3 

DIMENSION  HI  (12,12)  ,H  { 1 2, 12) ,5  { 12, 1 2)  ,G(12,  12) 
1  ,PHIT  (12,  12)  , 

*    Q(12,12),  G31  (120) ,  PKK(12,12) ,?KKH1  (12, 12) 
1  ,G11(120),G22{120) 

DIMENSION  IEEAD{10) , I WHITE (10) ,Y(120),YH(120) 
1  ,XH(120),G42(123) 

DIMENSION  DEL  (12, 12)  ,  A  (12r 12)  ,3(12,12)  ,D1  (12,12) 
1   ,02(12,12) 

DIMENSION  Z(12)  ,2(12)  ,GE  (12)  , XHP(120)  ,YHP  (120) 
1  , EY  (120)  , EX  (120) 

DIMENSION  DELI  (12,12)  ,PBI(12,12)  ,11  (120)  ,XKKM6  (1^) 


DIMENSION  ji&HI  (5)  ,  D  (12,  12)  ,  X  P  (2  5)  ,  Y?  (2  5)  rKJDLH3  (12  0)  , 
* V  ( 1 2, 1 2)  , f  (12,12)  ,  AI  (12,  12)  , 2KKH 1 ( 1 2) , TK ( 120) 
1  ,?KK22  (  I20)  ,P22  (120)  , 

*FT(12,  12)  ,  ET  (12,12)  ,XKK(12)  ,  XKXM 1  (12)  ,  ZXY  (12  J) 
1  ,E0B1  (120) 

DIMENSION  TIME  (12  0) , LAI (12  0)  ,LONG  (120)  ,ETIME  (120) 
1  , 3LAT  (120)  , 

*  E  L 0 &  G  ( 1 2  d )  ,  E VI N  D  ( 1 20 )  ,  X K M6  ( 1 20)  ,  X X M6  ( 1 20)  ,  PC N  ( 1 2 0) 
1  ,EOB2  (120) 

I  >l  I  EG S ^  IKM72  ( 1 00)  ,  T  K  .1 W  3  (100)  ,  I  K  3  2  U  (10  0)  , 

*  TKM12  (130)  ,irXo  (100)  ,IKP6  (100)  ,IKP12  (100) 
1  ,  IKP24  (103)  , 

*  TKP48  (100) ,TK?72  ( 100) 

J ATA    IQIT/1B$/,IY/' Y ' / ,1Z/ ' U' / , 1EE/1  t* / 

REAL    L  A I , I CN  G 

INTEGER    IYENO         /'Y'/ 
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I  N I E  G  E  5    I Y E N  1  /  •  Y  '  / 

NZ=    NO,     Of     OBSEEVED    VALUES  (SATELLITE; 
flZ^tlC-     Or     EESI     TRACK     VALUES 


liZ  =  50 
l4Z=5Q 

DC  132  1  =  1  , 4 
P  K  K  { I ,  I )  =  1  C  0  0 . 
132  PKKM1 <I,I}=1000. 
W  {1,1)  =  .000001 
.1  (2,  2)  =  .  CO 00 01 
2AE  OBSERVED  VALUES 

HEAD  (2,  11)  (TIME  {I)  ,  LAT  (I)  ,  LONG  (I)  ,PCN  (I)  ,I=1,NZ) 
11  FORMAT  (7X , F  6- 2 , F3 . 1 , 1 X , F4 -  1 , 1 X , F 1 . 0) 

FEAD  (2,  1  1)  (TIME  (I)  ,  LAI  (I)  ,  LONG  (I)  ,PCN  (I)  ,  1=1, NZ) 


:  READ  dEST  TRACK  VALUES 


READ  (3,  1  4)  (ETIME  (J)  ,  BLAT  (J)  , BLONG  (J)  , EWIND  (J) 

1   ,  j  =  1  #  a  z ) 

14  FORMAT  (6X,F4.Q,F4. 1, 1X,F4. 1,2X,F3.0) 

C  HEAD  (3,  11)  (ETIME  (J)  ,EIAT  (J)  ,BLOKG(J)  ,PCN  (J)  , J=1,MZ) 

C   11  FORMAT  (2X,4F10.2) 

C   11  FORMAT  (6X,F4.0,F4.1,  1X,F4.  1,2X,F3.0) 

C  ECHO  VALUES   ****REMOVE  NEXT  THREE  LINES  TO  ELIMINATE 

C       ECHO  PRINT  CHECK  *********** 

C 

WHITE  {8,373)   (TIME  (I)  , LA  I  (I)  , LON G  (I)  , E TIME  (T) 

1  ,3LAT  (I)  , 

*  3LONG  (I)  ,PCN  (I)  ,     1=1, NZ) 

C  WHITE  (8,  37  3)      (TIME  (I)  ,  LAT  (I)  ,  LOSG  (I)  ,  FCI^  (I)  ,1=1,  NZ) 
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J73  E08HAT  (2X,7F10.2) 
C      *P.ITE(9,373)   (TIHE(I)  ,LAI(I>,L0HG(I),ETIt1E(I) 

1  ,ELAI  (I)  , 
C     *      BLCNG  (I)  ,EWIND  (I)  ,  1=1, NZ) 
C      5? BITE  (9,374)  (Ilfl.E  (I)  ,LAT  (I)  ,IONG  (I)  ,E1IHE  (I) 
C      1  ,  SLAT  (I)  , 
C     *      PCS  (I) ,  I=1rNZ) 
C  J  7  4  FORMAT  ( 2  X  ,  6 1 1 0 . 2 ) 

c 


HIS  PHOGBAM  computes  the  following  kalian  filie:-. 


-1 


G(K)  =  E  (K/K-1)  *HI*  (H*P  (K/K-1)  *HT*£) 


?  (K/.<)  =  (I-G  (K)  *H)  *F  (K/K-1) 


P(K/K-1)  =  PHI*?  (K/K-1)  *PHIT+<2 


Q(I,J)  DEFINES  THE  COVARIANCE  OF  THE  PER  SAMPLE  RANDCfl 
EXCITATION  OF  THE  PROCESS 


R(I#J)  DEFINES  THE  BANDOJ3   (GAUSSIAN')  A  E  A  S  3  RE  .IE  NT  &')13l 
WHICH  IS  ADDED  TC  THE  CBSEBVAELE  SIGNALS 


Hi(I,J)  IS  THE  IDENTITY  MATRIX 


II  =  K  THE  DISCRETE  POINT  IN  TIME, THE  STAGE  of  THE  PRCCiiSS 
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C  PKK{I,J)     =    i }{K/K),      (-ov    EEF-OS    AT    K    GIVEN    K    SAMPLES) 

C 

c 

C  PKKK1(I,J)  =  P  (K/K-1)  ,  (COV  E8B08  AT  K  GIVEN  K- 1 

CSAtfELES) 

C   N  =  NUMBER  0?  ROWS,  H  =  NUMBER  CF  CCL.,CDS. 

C  ,  N  C  A  fci  D  .1  C  A  B 

C  Nj  =  MU3BEB  CF  ITERATIONS  CF  THE  FILTER 

C 

C 

q    *********  jj  g  y  NEWS  ************************************ 

fcRLTE (6,717 1) 
7171  FORMAT  (/  r  2X  #» *************  NEW  NEWS*** ****** *****'/, 
*2X,'THE  GTCF..1  I'BACK  INPUT  VALUES  ABE  A  71 AI  ABLE '/ , 
*-2X,  'THESE  ARE  ECHC  PBISTED  TC  THE  lEBHIKAL  AT  TEE'/, 
*2X, 'BEGINNING  CF  THE  PROGRAM  —  ALSO  SEE  LISTING'// 
*2Xf'IE  THE  PBCGRAfl  ENDS  NORMALLY  AM  INPUT  FILE 

1  aniVf 

*2X,'DE  PBGBUCED  THAT  dAY  EE  PRINTED  CF  USED 

1   FOB  INPUT'/* 

*2X,'USI8G    THE    SAME     EORilAI     STATEMENTS     10    READ    AS 

1  WESE'/, 

*2X,'USED    TC    WRITE    Cti    U  N  II    4.       5G    KITH    A    FEW 

1   JCDIr ICA1ICNS'// 

*2X,'THE  BRANCH  AROUND  THE  INPUT  CAN  BE  rJSED. 

1  PRESENTLY'/* 

*2X,'A  fILE  CF  THE  INPUI  DATA  IS  BEING  PRODUCED  AND 

1  BILL'/, 

*  2  X ,  '  3  E  FOUND  AS  — K  OUTPUT  A--  ON  YOUR  A  -  DISK-') 
C 

C   ****  IH3  FIRST  QUESTION  TC  THE  TERMINAL  --ASKS  IE 
C         AN  INPUT  FILE  IS  10  EE  USED 
C 
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h  5  IT  E  { 3 ,  2 7 ) 
27  format  1/2X,'  •/, 

*  t  ***-p***********  *********************************  I  / 

*2X,'D0  NCI  USE  ?A1  KEY  1C  EXIT  PROGRAM  UNLESS 

1  Y  C  J  H AMI'// 

*2X»T0    LOSE    YCfli:    INPUT    E  lui    THAT    IS    PRINTED    ON 

1    D8II    4.'/, 

*  i  ************************************************  »// 

*2X»FOB  CHECK  CUT  ACTIVATE  THE  GO  TO  833  STATEMENT 

1  c«v# 

*2X»LINE  135  TC  END  WITHOUT  ERROR1/, 

*i *************************************  *********** i  / 

*2X,'DO  YOU  KANT  TO  DEVELOP  THE  INPOI  FILE  [1/&)T) 

READ     ( 5  r  3  3 ) I fi  N  S 
33    FORMAT  (A4) 

IF  (IANS  .  NE.  IYEN1)  GO  TO  9  9  99 


7  7  7  7  F  C  R  M A  T { 1 H 1 ) 

WRITE  (8, 1234) 

1234  FORMAT ('ALL  INPUT  SHOULD  HE  FLOATING  POINT  E<lEPT  ') 
WBIIE{9,1235) 

1235  FORMAT  ('  WHERE  OTHERWISE  SPECIFIED  ',/) 
10  1b  FORMAT  (1H0,5X,5H  K  =  ,  13  ,/o  '< ,  5HG  A I  N5) 
19  3  2  W  B I T  E  ( 3 ,  1 0 ) 

10  FORMAT  {/,5X, '  THE  DISCRETE  KALMAN  FILTEF1) 

WRITE(8,30) 
30  FORMAT (5X,  •      ENTEF  THE  O^DES  OF  THE  SYSTEM 
1  (UP  TO  3)  .  ') 

S  £  a  D  ( 5  ,  4  0 )  N 
4  0  FORMAT  <I1) 
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ND=  1  2 

a  d = 1 2 

DC  7838  1=1,  N 

Z  (I)=0. 

DC  7  8  9  3  J  =  1 , N 

HI  (I,J) =0. 

HI  (I,I)  =  1. 

? HI  (I, J) =0, 

PHI  (I,I)=1. 

H (I,J)=0. 

E  { I ,  J )  =  0 
7898  A  (I, J) =0. 

WBITE  (8,4  10) 

^RITE  (8,7771) 
777  1  FOBHAT  (5X,  »CO  YO'J  WANT  1C  COHPOIE 
1     LINE  F?C.:V  , 
+  1  X  ,  •  A  S  3  ?  ' ) 

aEAD(5,7772)IAN 
7  7  7  2  F  C  B  ;1  A  T  { *\  J ) 

IF  {IAN  .EQ. IPH)  GO  TO  7773 

IF  (IAN  . NE.IPH) GO  10  15 

7773  WSITE{3,7774) 

7774  FORMAT  (//,5X  ,' ENTE?  THE  A  «ATEIX») 
DO  7775  1=  1, N 

DC  7  77  5  J=1,N 
WHITE  (3,777c) I, J 

77  76    FORMAT  {5X,  •  A  {»  ,H#  ','  /H/M  =') 

READ<5,7778)A(I,J) 
7  77  8    frOBMAI{F10.3) 
777  5    CO NT  IN OS 

WHITE  (9,77  79) 

WBITE(8,7779) 
77  79    FCP.-1AT  (5X,  '       THE    A     .1 A  IRIX  •) 

DO    7 78  J    1=1  ,N 
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BPIT2{*,90)  (A  {I, J)  ,J=1,N) 
77-0  WRITE  [8,90)  (A  (I, J)  ,J=1,N) 
C      WRITE(8,4422) 

C4*z2  FORMAT  (5X,«   ENTE5  THE  DIMENSION   CF   A  ') 

C 

C      Vi  IS  CGN3TEAINED  1C  (1,1)  HERE  6  2  KK  Q455GS  0270  J 
C 
C 

C      F5 A 0(5,4425)1 
0^423  FORMAT  (11) 

1=1 

WRITE  (8,7781) 
77  5  1  FORMAT  (5X,  '  ENTER   -HE   E  MAI  SIX   •) 

DC  77-32  1=1,8 

DO  7732  J=1,L 

WRIT £{8,7783)1, J 

7733  FORMAT  (5X, fE<« ,I1,», f ,11,*)=') 
SEAD(5,7778)B(I,J) 

77  3  2  CONTINUE 

w'SITF.  (9,77  84) 
53 HIT"5!  (8,77  84) 

7734  FORMAT  (5X, '   THE  0  MATRIX    ») 
DC  7735  1=1, N 

W3ITE<9,90)  (B<I,J)  ,0=1,1) 
770  5  fcRlTE  (3,  30)  (E  (I, J)  ,  J=1  ,L) 

GO  IC  90  10 
15  w'3ITE(3,5G) 

50  FORMAT  (5X# *  ENTER  IRE  ELEMENTS  OF  IRE  I  RAMS  I  ION 
1  S  A  T  RI X  —  P  B I » ) 

DO  1  1=1, N 

DO  1  J=1,N 

BRITE  (8, 60) I, J 
b 0  F0RMAT(5X,*PHI 

READ  (5,70)  PHI  (I, J) 
7u  FORMAT  (F  10.0) 


,11,', »,I1,»)= 
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1  CONTINUE 
«RITE(8,30) 

BO  FORMAT  {'  3'  ,  5X,  •  THE  PHI  BATEIX  (TRANSITION  MAIBIX)*) 
DC  2  I = 1 , N 

2  SRITE{8,90)  (PHI  (I, J) , J=1 ,N) 
90  FORMAT  ( 1 P 7 £  1 1-4) 

9  1  FORMAT  <1P7Z1 1,U) 

£*##*****#*#$£****•*#**£******** 

c 

C      GO  TO  883 

10  00  WPI'IE  (dr  10G) 
100  FORMAT  (5X,  'DO  YOU  WANT  TO  CHANGE  ANY  ELEMENT  OF  IHE 

1  HAIEIX?1) 
900  8EAD{5,110)IAN 
110  FORMAT  {A3) 

IF  (lAN.SQ.IZ)GOTO  19 
IF  (IAN.NE.  IY) GOTO  1000 
SPITE  (8,12  0) 
120  EORMAT  (5X, 'RH.ICH  ELEMENT  OF  IHE  MATRIX  00  YOU  KANT 
1  TO  CHANGE?' ,/, 

+5X, 'ENTER  AS  IJ;  WHERE  I  IS  IHE  ROti  AND  0  IS  IHE 
1  COLUMN- •) 
READ  (5, 130)  I, J 
130  FORMAT  (211) 

a  R I T  E  { 8  ,  6  0 )  I ,  J 
BEAD  (5,70)  PEI  (I, J) 
aRITE(8,30) 
DO  3  1=1, N 

3  MRITE  (8,90)  (PHI  (I, J)  ,J=1,N) 
SEIIE{8,14Q) 

140  FORMAT  (5X, 'ANY    OTHER    CHANGES?') 
GCTO    900 

19    WRITE  (9,  14  1) 

141  FORMAT  (' 1« ,5X, 'DISCRETE  TIME', 
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MX^KAI'IAN     FIITES     PPGGBAS') 

51  BIT  £  ( 9  ,  1 4  2)  (N  AME  ( Ij  ,1=1,5) 
142    FCRJSAT  (6X, 'FEOBLEH    IDENTIFICATION  :  *  ,  5X  ,  5  A4) 

WBITE(9,143) 
1  4  3    FGRdAT  ( *  0  »  ,70  { '  *« )  ) 

WRITE  (9,30) 

DC     1111     1=1, N 
1111     WRITE {9 , 3 1 )   ( 5 H I  ( I , J)  ,  J=  1  ,  N ) 

WRITE  (3,  4 1G) 

HRITE{8,150) 
150    FORMAT  (5X, 'ENTER    1H£    DIMENSION    CF       BANDOM    INPUT 
1     V2CTQB     (8)  ') 

FEaD  (5,40) I 

SRITE  (8,  160) 
1o0  FORMAT  (5X,  »ENTER  THE  ELEMENTS  OF  THE  DISTRIBUTION1, 
+  1X,  '  MATRIX— GAMMA.  '  ) 

DO  4  I =  1  ,  M 

DG  4  J=1 , I 

8BITE(8,170)I,J 
170  FOBMA1  (5X,  » GAS  HA  ( '  ,  I  1  ,  '  ,  '  ,  I  1  ,  » )  =  ') 

5EAD(5,70)DEL(I,J) 

4  CONTINUE 
WRIT E{8, 180) 

160  FORMAT  ('(J*  ,5X,  'THE  GAMMA  3aT?.IY  (DISTRIBUTION 
1  H  A I R I X )  » ) 
DC  5  1=1, N 

5  it  R  IT  E  { 8 ,  9  0 )  (DEL  (I , J)  , J  = 1  , L ) 
110  0  »aiTE{8,100) 

9  10  BEAD(5,110)IAN 

IF  (IAM.EQ.  IZ) GOTO  29 
IF (IAN. N 2. 1Y) GOTO  1100 
WRITE  (8, 12  0) 
EEAD(5,130)I,J 
WRITE (8, 17 C) I, J 
BEAD(5,70)GEL(I,J) 
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w3ITE  (8,  18  3) 
DC  6  I = 1 , N 
6  KRITE(8,90)  ( DEI  (I , J)  , J=  1 ,1 ) 
WRITS  ( 3 , 1  4  C ) 
GOTO  9  ]0 
2  9  SfiIT£(9,180) 
DO  1112  1*1, a 
1112  WRITE  (9,91)  (DEI {I, J)  , J=1 ,1) 

»5ITE(8,41Q) 
90  10  WRITE  (3,  393) 
330    F0R3AT  (5X,  •ENTER  THE  ELEMENTS  OF  THE  COV  OF  B») 
C 
C  NOTE  %       IS  CONSTRAINED  10  (1,1)  3  KK  17  10  5  '4450 

W  (  1 ,  1 )  =  .  J  0  3  G  0  1 
C 

7c  DO  21  1=1,1 
DO  21  J=1,L 
WRITE  (8, J 10)  I, J 
3  10  FORMAT  (5X,  '«  ( '  , I  1 ,  ' ,  ■  ,  1 1  ,  ' )  =  ■  ) 
REAL)  (5,70)  'a  (1,0) 

21  CONTINUE 

7  7  WRITE  (8,3  21) 
WHITE  (9, 32  1) 
321  FORMAT  (»()•  ,5X,  'THE  COV  OF  W») 
78  DO  22  1=1,1 

22  WRITE(3,90)  (W(I,J)  ,J=1,  L) 
1u00  WBITE{8,100) 

94  0  ?EAD(5,110)IAN 

IF  (IAH.EQ.  11)  GOTO  69 

IF  (IAN.NE. IY) GOTH  140  0 

«EIT£  (8,  120) 

READ  (5,  130) LfJ 

WEITE  (8,310) 

R  E  A  D  ( 5  ,  7  0 )  ii  (1,0) 

M  R I T  E  ( 8 , 3  2  1 ) 
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23 

89 

1115 

190 
195  F 


20  0 


CO  2  3  1=1, L 

WRITE  (8,90)  (*  (T,  J)  ,0=1,1) 
WBITE(8,140) 
G U  T  0  9^0 
DC  1115  1=1,1 
WRITE  (9,91)  (W {I, J)  ,0=  1  ,1) 
WBITE(8/190) 

FORMAT  (5X,  'ESI ER  'IRE  CEDES  OF  H,  I.E.  3=?») 
li  E  A  D  (5,40)  B 
WBITE(8,195) 

FORMAT <5X, 'ENTER  IEE  EIESENTS  OF  THE  OBSERVATION 
1  MATRIX — B.») 

co  7  i  =  1 ,  ;i 

CO  7  0=1,11 

WHITE  (3,200)  1,0 

F  C  B  3  AT  ( 5  X  ,  '  fl  (  •  ,  1  1  ,  •  ,  '  ,  I  1  ,  »  )  =  '  ) 

READ  (5,70)  H  (I, J) 

CONTINUE 


aSITE  (8,210) 
210  FORMAT  (,0»,5X,  «THE  B  HAS  BIX  (OBSERVATION  3ATRIX)*) 
C  C  8  I  =  1  ,  A 
3  BBITE  (8,90)  (H  (I, J)  ,J=1,N) 
1z00  WRITE  (8,10  0) 
920  BEAD  (5,  1  10)  IAN 

IE  (IAN. EQ.IZ) GOTO  39 
IF  (IAN. NE. I Y) GOTO  1200 

.«  r  I  I  J.  (  0  ,  U  C  ) 

HEAD  (5,  1  JO) 1,0 
WRITE  (3,200)  1,0 
BEAD(5,70)B(I,J) 
WBITE(3,210) 
CO  9  1=1, M 
9  WRITE(8,90)  (II  (1,0)  ,J=1,N) 
WBITE(8,140) 
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GOTO  9  20 
J 9  H SITE  (9,210) 
DC  1113  1=1, fl 
1113  »  5  IT  I  ( 9 , 9  1 )  ( B  { I , J)  , J  =  1  , fi) 
WEITS{9,143) 
WHITE  (8 , 270) 
2  70  FGBflAT  (5X,  'E81ER  THE  EIEflENTS  Or  THE  3EASUBEMENI 
1  NOISE', 
+1*, 'COVABIANCE  MAIHIX — F   •) 

12  DC  116  1=1, a 
DC  116  J=1,M 
WHITE  {8,28  0)1, J 

20  0  F0RflAT(5X,  »fi { ' ,1 1 , ' , » , 1 1 , ' ) = * ) 

BEAD  (5,70)  B{I, J) 
116  CONTINUE 

WHITS  (3,290) 
WHIIE{9,290) 
293  FOfiflAT  ('0  '  ,5X,  ■  IHE  H  BATBIX  (HEASDHEHENT  •  , 
+  1X/'NCISE  CO VARIANCE  rt  A I ^ I X )  • ) 
7  4  DC  17  1=1,  K 

17  WBITE(3,90)  (F  (I, J)  ,  J=  1  ,  M) 
130  0  HBITE(8,  100) 
930  HEAD  (5,110)  IAS 

IF  (I AN. EC;.  IZ)  GOTO  7  9 

IF (IAN. NE. I Y) GOTO  130  0 

WRIIE{8,120) 

BE  A  D(5, 130)1, J 

»BITE(8,230)I,J 

READ  (5,70) B (I, J) 

»BITE(3,29C) 

DC  1d  1=1, a 

13  WHITE  (8 , 90)  (H  (I , J)  , J= 1 , M) 
aBITE  (3,  140) 

GOTO  9  30 
79  DC  1114  1=1, fl 
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11  14  KRIIE(9,91)  {?  (1,0)  ,J=1  ,S) 
WHITE(8,410) 
H  B 1 1 E  ( 8  ,  J  3  0) 


5  0  FOR MAI  (//,5Z,  •  ENIE 


FK SM1  (1/0) 


) 


4/  CL  J4  1=1,  N 
DC  34  J=1,N 
8RITE  (8,340)  I, J 
3 4  0  FOHHAT  (5X, "EKKH1  ( ■  ,  I  1 , ' , • , 1 1 ,  ' )  =  •) 
READ  (5,70)  PKKM  1  (I,  J) 
34  CONTINUE 
4  3  WHITE  (8 #  35  1) 
m  8 1 T  E  (9,351) 
3  51  F  C  3  A  A  T  ( '  0  •  ,  5  X  ,  '  ?  rC  K  a  1  (  1  /  0  )         '  ) 
4  4  DC  3  5  1=1,  N 

3  5  H BIT E  (3 , 9 0 )  (PK Kfi 1  (I , J)  , J= 1 , N) 
150  0  SHITE  (8, 100) 
950  F.  £  A  D  {  5  ,  110)1  AN 

If  (IAN.Ee. IZ)  30TO  51 
IF  (IAN.N3. IY) GOTO  1500 
WFITE(8,12G) 
READ  (5,130)1,0 
WBIT£(8,340) I ,  J 
.READ  (5,70)  PKKH1  (1,0) 
SHITE  (8, 351) 
DO  3d  1=1, N 
3  b  U SIT E  ( 3 ,  9 0)  ( PK Xii  1  (I  ,  J)  ,  J=  1  , N ) 
HBIIE(8,140) 
GOTO  950 
51  DC  1116  1=1, N 
1116  WHITE  (9,91)  (EXK.11  (I,  J)  ,0  =  1  ,  U) 
WRIIE(9,143) 
KEIT£(8,246) 
2^o  FOPH&T  (5X,  •EN1EE  TH5  80MEEB  OF  IHE  POINTS  TC  3E 
1  PERFORMED. •  , 
Vf8X#,(<100)  SHIS  IS  AN  INTEGEH  VALUE') 
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3  HAD  (5,24  7)  NM 
2  4  7    F  C  R X  A  T  { 1 2 ) 

IT=  (TIME  (1)- 1400)  *  (.01) 
IS=  (TIME (2) -1400) * (.01) 
1-24*11 

S=24*IS 

T=IIME  (1)-IT*100-1400+T 

S=1IHE  |2)-IS*10Q-1400  +  S 

DT=S-T 

27^3  CONTINUE 
JK=1 
HRIT£(8,7865)JK 

7365  FORMAT  (14) 

i«RITE{8,51  1)  N,M,L,  ND,MD,LD,NN,DT 
511  FORMAT  <2X,2f!N=,  15,  5X,2HM=,  15,  5X,2HL  =  ,  15,5 X/3HND= 
1  ,I5,5X,3HMD,I5# 
*5X,  3;lLE=, 

&I5,5X#3HNN=,I5,5X,3HEI=,F10.4) 
WRITS  (3,533) 
533   F0RMA1(/'    MATRIX  R  »  ) 


3017 


544 


3013 


555 


019 


DO  30  17  1=1, a 

*RITE(8,90)   (R(I,J) ,0=1,8) 

SRITI  (8,54  4) 

FORMAT  (/'    MAT R IX  Q  • ) 

DC  3  0  18  1=1, N 

HEITE{8,90)  (Q(IfJ) ,0=1, *) 

» RITE (3,5  5  5) 

FORMAT  (/'    MATRIX  P K K 31  ' ) 

DO  3019  1=1,  SI 

WRITE  (3, 90)  (PKKM1  (I, J)  ,J=1, N) 

EF  (IANS.NF.IFH) GO  TO  6789 
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11=  (TIME 
15=  (TIME 
1=24*11 


( 1 )  - 1  4  0  0 )  *  (-01) 
(2J-1400)  *  (.01) 


S  =  2  4  *  I S 

I  =  T I  -A  E  (1 )  -IT*  1 00-  1 400+1 

S=TIME  (2)-IS*10  0-1400+S 

DT=S-1 

Ja  =  2 

WRITE(8,7365)OK 
9753  CONTINUE 

WHITE  (8,78  78) Dl 
78  78  FORMAT  (5X,F1 0.3) 

CALL  PHIDEL (BT, N,L, A, E,PKI, DEL, D1,D2,N D,MD,LD; 

WRITE  (8,666) 

FORMAT  {/'      PHI       • ) 

DO  3020  1=1, 8 

WRITE  (3,9  0)  {FHI  (I,J)  ,0=1, N) 


6739 
6oo 

3020 


mP.IT  E  (8,777) 
77  7    FORMAT  {/ «      DEI        • ) 
DO  3  0211=1  ,N 
30.:1    WRITE  (8,90)  (DEL  (I, 1) ,0=1 , L) 

CALI  TRANS  (BEL, N , 1 , DELI, ND, HD) 

CALI  PROD  (DEL, DELT,N,  1,8, Q  ,8E,MD,LD) 

CALL  CONST  (ti  (1,  1)  ,Q,N,N,Q,ND,  ID) 

BSITE  (3,54  4) 

DO  3  02  5  1=1,  N 

3025  WRITE  (3,90)  (Q  (1,0)  ,0=  1,8) 
3RIIE(3,iJ44) 

444    FORMAT  (/  '  II      •  ) 

DO  J026  1=1,  .1 

3026  WRITE(3,90)  (3  (1,0)  ,0=1 ,N) 
JK=3 

MRITE  (8, 7865) JK 
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iifiITE(<*,5111)K,!!,L,HD,?1D,LD,NN.,DT 

5  111  E  G  E  A  A I  {  7  I  4  ,  :  1  C  .  4 ) 

DC  5  32  7  1=1, N 
3327  KfiITE(4,90)   (A  (I,J)  ,J=1,N) 

DO  5  3  28  1=1, N 
3323  isSITE(4,90)  E(I,1) 

DO  5329  1=1, N 

5329  »BITK(4f90)  (0.(1, J)  ,J=1,N) 

DO  523C  1=1, S 

5330  5?BITE(4,90)   (EKKtf  1  (I,  J)  ,  J=1  ,  N) 
DC  6327  I =  1 , fl 

o3z7  hFITF  (4,93)   (fl  (I, J)  ,J=  1  ,  J) 
DO  6  927  1=1,8 

6  9  2  7  S  BIT  E (4 , 90 )   ( KI  (I , J)  , J  =  1 , N ) 

DO  6  3  26  1=1, H 
6326  »[jiTE  (4,90)   (E  ( I ,  J)  ,  J=  1  ,  .1) 

rwRITZ  (3,777  7) 
9999    CONTINUE 

IF  (IANS.EQ-IXEN1) GC  TO  7234 

WHITE  (8, 1928)  IANS 
1926  7C?:!AT(5\,'  XXXXXXXXXXXXXX  •  ,1W) 

JK=3o 

WFITE{3,7865)JK 

READ(4,5111)N,M,L,ND,MD,LD,NN,D1 

JK=  J  9 

» BITE  (8, 7865) JK 

n  BIT E  (8,  5 1 1 1 )  N  ,  3 ,  L,  SD ,  a C  ,  ID ,  N  HI ,  CI 

JK  =  4 

WHITE  (S,76c5) JK 

DC  723j  I =  1  , H 
7285  HEAD(4,91)   (A  (I ,  J)  ,  J=1  ,  N) 

JK=4 32  1 

W5ITE{8,7865)JK 
DO  7  28  6  1=1, N 
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7  2  j  c  5  E  a  0  { '4  ,  9  1 )  E  ( 1 ,  1 ) 

DC  7  z  3  3  I  =  1  ,  N 
7z3d  READ(4,91)   (Q  (I,  J)  ,J=1,N) 

DC  7239  1=1, N 

READ  (4  ,  9  1)   ( F  K  K  S  1  ( I  ,  J)  ,  J=1  ,  N) 
7239  WRITS  {3, 90)   (PKKH1  (I ,J)  , o=  1  , N) 

DC  73  2  6  1=1, H 
732C  PZAD(4,91)   (H  (I,J)  ,J=1,N) 

DO  7j29  1=1, N 
7329  SEAD(4,91)   (HI  (I  ,  J)  ,  J=1  ,  N) 

DO  73  2  7  1=1, a 
732  7  HEAD  (4,91)   (B  (I,  J)  ,  J=  1  ,  M) 

IT=  (TIME  (1)-150Q)  *  (.01) 

13=  (TI/!'7  (2)  -  1500)  *  (.01) 

1=24*11 

S=24*I3 

i= II. IE  (1) -11*10  0-1500  +  1 
S  =  II:1£  (2)  -IS*  100-  1500 +  S 
D I  =  S  -  i 
92  5  3  CONTINUE 

R  ( 1 ,  1 )  =  -  0  0  0  0  0  1 

JK  =  D 

WRITE  (d, 7dfc5) JK 

CALL  TRANS  (  D  E  L  ,  N ,  1  ,  0  E  L 1  ,  N  D  ,  _kl  D ) 

CALL  PROD  (DEL, DELT,N, 1,N,Q  ,ND,aD,LD) 

CALL  CCNS'I  (w  (  1  ,  1)  ,Q,N,N,Q,ND,8D) 

CALL     PHID2I  (CI  ,  N,L,  A,  E,  Fill,  DEL,  D1  ,  02,  ND,  '10,  ID) 
72  34    CONTINUE 

D  D I  =  D I 

a  R  1 1 E  ( 3  ,  9  1  5  z ) 
9162  FORMAT  (5X,  '  YCU  All  AT  DO  2222  •   ) 
C  DLDLDLDLLDDLDDLDLDLDLDLDDLDLLCLDLDLDLDLEDLJ 

EB=  150  0  +  11* 100 

KKK=  1 
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DO  2  22  2   K  =  1,  v''i 


C       ***  CALLS  SUBROUTINE  TO  CALCULATE  JULIAN  I L A E 
c       ***  FCn  EVEEX  ST08M  POSIIICN  AND  E7IEY  fi  HOU8S 

ITISE=INI  (TIME  (K)  ) 

IDAY=III JE/100 

IHOUB=ITIM£-  (IDAZ*100) 

IF  CIHOOe.EC-0) IHOUE=24 

3CD6  =  :iOD  (IHCUB  ,6) 

CALL    JOTIfilE     (ITIME/JOLBEJ 
C  IF  (MOD6.EQ-0) GOTO     1986 

HTIHE=IN1  (EC) 

MDAY^TI  -IE/100 

N  D  A  Y  =  11  D  A  Y  *  1  0  0 

£!HOlJH=aTiaf-MEAI 

IE  (HHOUR.  31.  24)  ilTIH£=ITIHE 

CALL    JUTIME     lflTIHEfHJUlflB) 
193b       CONTINUE 

»SIIE{8,  19  89)  JULEIE,III3E 
1989       FGI-MAT  (///,'     JULIAN     HOUB    IS     ',19,'     , ACTUAL    IId£    IS 
1        ',1b) 

SHITE  (3,  19-84)  HT:iaE,HJaifl£ 
1984       F05aAT{«     30DU1A    6    TI;1E  = 

1  JULIAN  H B  = ' , 1 9 ) 
198  5   CGNTINUE 

C       ***  END  JULIAN  TIME  hCUTINE 
C 
C 
,1RB 


•  , 15,  '  rCOBBESPONDING 


***  CALCULATE  HODULA  6   JULIAN  II JS  FOE  1KB, IKH, XK P 


C 

C 

c 


***  NODULA  o  EOT  JULIAN  IIM 
IF  (MHOUR.EC- 18) GOTO  3187 

IKH72  (K)=HJ0LHB-7  2 

IKM48  (K)  =.'JJULHR-43 
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T  K  *  2  4  (  K )  =  :•]  J  rJ  I  H  7  -  2  a 

T012  (K)  =HJDLHE-12 

T  K  M  6  ( K )  =  3  J  D I H  S  -  6 

TKP6  (K)  =1J'JlcE.  +6 

TKP12(K)=MJCJIHR+12 

IKP24  (K)  =3JUIHR+24 

TKF43  (K) =MJUIHR+48 

lK27z  (K) =aJ0Lfifi  +  72 
3167   CONTINUE 

C      ***  END  JULIAS  IlflE  BOUTIN E 
C 

IF  (PCS  (K) -5. NE.O)  GO  TO  3133 
3  {1,  1)=.25 

R  (2,  2)=. 25 
GO  TO  3134 
31j3  IF  (PCN  (K) -3- NE.O)  GO  TO  3135 
3  (1,1)  =.0625 
R (2, 2) =.0  625 
GC  IC  3134 
3135  IF  (PCN  (K) -2. NE.O)  GG  TO  3134 
B{1,1)=.0312 
R  (2,2)  =-0312 
3  134  CONTINUE 

SHITE  (8,544  5)  IlJiE(K) 
544  5    FORMAT  (////// 5  OX , 5  HI IflE  = , r 10. 4) 
SRITI  (4,54  45)  I  IKS  (K) 
WRII£(8,  *313)K,3E,B  (1,1)  ,C  (1,  1) 
9  313  FO  2  N  AT  (  3  X  ,  •  K=  '  ,  1 3  ,  5  X  ,  '  b  5  =  «  ,  ?  8  .  2  ,  5  X  ,  •  B  =  •  ,  F  7  .  4  ,  5  X  , 
1   »Q  (1,1)  =  ' F  10.  4) 
*PITE(4,9313)K,23,5  (1,1)  ,Q  (1,  1) 
*  RUE  (8,3  1  2)     PCN  (F)  ,DT,8  (1,  1) 
3  13  FORMAT  (3X,  '  PCS  (K)  =  •  ,F10.  4,  3X,  «  LT=  »  ,F6.  2,  -10X,  *  B  { 1 ,  1) 
1  =  '  ,  F  1  'J .  4 ) 
'A  RI T E  {  4  ,  3  1  2)  PC  N  (K)  ,  DT  ,  w  (1,1) 
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3  1 2  9  Q  ( 


DO  3129  I  =  1 ,  X 
DC  3  129  J=1,N 

Q(l,3)     = 


U  (J.,  J)  -  u. 

5(1,1)  =  (ET**4, 

Q  (2,2)  =  (E 1**4/4)  * 

-  12,  3 )  =  n  t  *  *  ?  *  i 


Q(4,4)   = 


^  \  -r  ,  -r  /         -        u  J. 

WRITS  (3,799) 
7  99  F0B3AT  (/  * 

DC  3123  1=1,1- 
3123  WRITE  (8,90)  (C 


/4)  * 

,/4)  * 

Dl**2*W(2r2) 
2*W(2,2) 


(1,1) 
a  (2,2) 


MATRIX   Q      ' ) 


DC  3123  1=1,  N 
WRITE(8,90)  (C  (I, J)  ,J=1,N) 
W (1, 1) =.00  0001 


031 
9  (2,2) =.000001 

CALL  PHI DEL  (ET,M,I 
dRITE  (8,97  9 


, A, E#PHI, DEL 


D1,D2,ND,3D,LD) 


2) 


9  79    FOR 


Pril    •) 


■  r  -■  / 

J=1,N) 

i,  g, a, phi, a, n,s,g, hi, hd,md 


MAI  (/  • 
WRITE  (4,97  9) 
DC  3579  1=1, H 
wRITE{4,90)  (PHI  (I, J)  ,J=1,N) 

579  WRITE  (8,9  1)  (Efil  {I,  J)  ,  J=  1  ,  N) 
CALL  GAIN  (PKK,PKKM1,Q,B,P 
1      ,LD,K) 

WRITE  (4,656) 

WR I TE  (3  ,  6 5 ' 

FORilAT  (/  ' 

DC  30  2  3  1=1, X 

W  P  IT  E  ( 3  ,  9  0 )  ( PKK  ( I , J )  , J  =  1  ,  X ) 

WRITE  (4 , SO)  ( PKK  (I  ,  J)  ,  J  = 1  ,  X ) 

CALL  PR-CD  (PIiI,XKK,N,ft,  1 ,  XX  K  M  1  ,  XD  ,  ED  ,  ID) 

CALL  PRCD  (H,  \KKM1,  X  ,  >1  ,  1  ,  2  K  K  .1 1  ,  XD,  MB,  3D) 

WRIT E(3, 8810) 

P0R3  AT  (/  •     ZKKM1       M 


56) 


65o 


3023 


PKK 


38  10 


—  —  —  \-#  —  —  '/ 

F0R3  AT  (/  •     ZKKM1       * ) 
WR ITE  (8 , 90 )  (Z KKfl  1  ( J )  , J= 1 , A ) 

WRITE  (4,  90)  (ZKKfll  (J)  ,0  =  1,11) 

WRITE  (d,  38  19) 
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S819  FOa.iaT(/»     LAT  (S)  ,LOSG  (K)      '  ) 

C  WHITE(3,90)IAT(K),IGNG(K) 

Z(1)=1A1  (K) 

Z  (2)  =LOHG  (K) 

W  H  IT  E  { 4  ,  8  8  1 1 ) 

WRITER, 30)  (Z  (J)  ,J=1,M) 

WHITE  {3 ,88  11) 

WRITE  (8,90)  (Z(J)  ,J=1,«) 

CALL    5U3  (Z,ZKK21  1  ,  2!  ,  1  ,  E,  8  D,  HD) 

WRITE{8,5445)     T I  fl  £  ( K ) 
C  KRITE(3,S810) 

WBITE(3,3029) 
3029  F0  3M5T  (/     '  E  *^******^**^*j(c************* 

WRITE  (8,90)  {  £  (J)  ,J=1,N) 
IF  (K.LE. 1)  GO  10  2  20  4 
GME=  (PKKtfl  (2,2)  +5(1,1))  **.5 
IF  (AES  (E  (2) ) -GATE  .11.0.)  GO  10  2203 
G  (2,2)  =0.5*  (1. 2  +  G  (2,2)  ) 
W  (2,2)  =10000.  *W  (2,2) 
G  (4,2)  =0.5*  (0.  333  +  G  (4,2)  ) 
C      PKKH  1  (  1 , 1 ) -2 *P  KKfl 1  { 1 , 1 ) 
C      PKKMl(2,2)=2*PKKi3l<2,2) 
C      PKKB1  (3,3) =2*PKKM1  (3,3) 
C       PKK/i  1(4,4)  =2*PKKH  1  (4,4) 
C      WRITE  (4,9  192) GAIL,  2  (2) 
H 3 IT E  ( 8, 9 192) GATE, E  (2) 
2203  CONTINUE 

G AT E= (PKKfl 1  ( 1 , 1) +B  ( 1 ,  1 ) ) ** . 5 

IF  (ABS  (£  (1) ) -GATE  .LT.O.)  GC  IC  2204 
G(1,  1)=0-5*(1.2+G(1,1)  ) 
W(1,1)=1v)000.*fc(1,1) 
G{3,1)=0.5*(0.333  +  G(3,1)) 
C      ?KK£31{1,1)-2*PKKM1<1,1) 
C      PKK.11(2,2)=2*?KKH1(2,2) 


') 
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,F10.4,9X,»E  (2)  - 


C      PKK31  (3,  .3)  =2*PKKH  1  (J,  3) 
C      PKXd  1  (4  ,  4)  =2*PKKi1 1  (4 ,  4) 
C      WPITE  (4,9191) GATE, E (1) 
WRITE  (J,  9  19  1)  GATE  ,E  (1) 
9191  F0F.3&T  (9X,  'ERFCR  JT  GATE.  GATE=  »  ,  F  10.  4,  9X ,  •  E  ( 1 )  = 

1   ,  F 10-4, 'XXX') 
9  192  FCR.1AT  (9X,  'ERROR  GT  GATE.  GATS- 
1  , F10.4, 'XXX') 
C 
C 
C 

c 
c 
c 
c 
c 


G1  1  (K)  =G  <  1  ,  1) 

331  (K)=G   (3,1) 
DC  3022  1=1, U 

WHITE  (6, 90]  (PKKH1  (1,3)  ,J-1 ,N) 
30  22  wRITE  (4, 90)  (PKKS 1  (I , J)  , J= 1 , N) 
2204    WRITE{4,99) 

WHITE  (3, 59) 
99    FORMAT  (/  '    HATHIX  G   •) 

DO  3024  1=1, N 

WRITE  {8 r 9 G)  ( G  (I ,  J)  ,  J=  1 , fl ) 
30  2  4    WRITE (4,90)  (G(I,J) ,  J=1  ,fl) 

WRITE  (3,  90)  (ZKKfS  1  (J)  ,J=1  ,3) 

H RITE  (3,88  11) 

WRITE  (4,38  11) 
63  1  1    FORMAT  (/  •      Z      •) 

*3IT£  (3,  90)  (Z  (J)  ,  J=  1,i1) 

*RIIE(4,90)  (Z  (J)  ,J=1,H) 

CALL  PROD  (G,  E,  .<f  M,  1  ,  G  Z  ,  iN  D  ,  :1D,  LD) 

WRITE  (4, 90)  (GE  (J)  , J=1,N) 

CALL  ADD  (XKXH1,GE,K,  1  ,XKK,ND,3D) 
:        CALL  PRCD  (PHI,XKK,N,N,  1 ,  XKKfl T, ND,3D , LD) 

rfRITE(8,8011) 

W  ?  I T  E  (  4  ,  o  0  1  1  ) 
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3  C  1  1     F  C  ''  3  A  7  (/  '     X  K  K  »  ) 

SSITE  (8,  90 )  (  X  7  E  ( J )  ,  J=  1 ,  N ) 
«2ITf  (4#  30)  (XKK  ( J)  , J=  1  , N) 
311{K)=G  (1,1) 
GJ1  [K)=G  (3,1) 
C        ?1  1  (K) =PKKM1  (1,1) 
C        ?K1 1  (K) =  EKK  (1,1) 
Y  H  (  k  )  =  X  K  K  (  l ) 
X  H  (  K )  =  X  K  K  (  2 ) 
KTIHE=IN1  (TIME  (K) ) 
RKTIMS=FICAI  (KTI3E) 
x?  (HKTI:1£.  HE.TI3E  (X)  )  GOTO  83  13 
KJOLHB  (K)=J0LHJ 
L         EOB1  (K)  =YH(K)-LAT  (K) 
C        EOB2  (K)=XH  (K)-LONG  (K) 

36  13    CONTINUE 
C        WRITI{4,8812) 
C        WHITE  (3, 88  12) 
66 12  FORMAT"  {/ i        XKKM1     '  ) 

fcBITE  (3,  '30)  (XKKfl  1  (J)  ,J=  1  ,N) 
S  HIT  E  { 4  ,  9  0  )  (XK Kli  1  ( J)  ,  J= 1 , N ) 
YHP (K) =XKK«1  (1) 
XHP  (K) =XKK«  1  (2) 
8 3=1500+  11*100 
96  35  DC  968  4  1=1,4 

3686  IF (TIME  (K) -EB-LE.0)  GC  IC  9937 
Et=33+6 

IF     (BB-TIflE  (K)  .LI.  0)  GC  IC  36oo 
9997  IF  (33-TI'1E  (K)  .GT.  0)  GO  TC  9o99 
Y  K  3  b  ( X  *  K )  =  X  K  K  (  1 ) 


C 


XXM6  (KKK)  =XKK  (2) 
GO  TO  9  693 
9  6  9  9  CC=BB-  ( 1 500  +  IT *  1 00 ) 
WHITE (4,56  56) CC,bB 
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WRITE  (8,563c) CC, 3E 


5  c  5  o  F  G  R  cv 


A I  { 5  X ,  •  i 

13  (CC-2  4. EC 
IF  (1-4. LI. G 


:  c  = » ,  f  1 

))  GO  I 


n  u  » 


153='  ,F1J.  4,  ' 


30K  '  ) 


9  7  0  0 
re  9o94 


i  ir  ( i  -  4  .  L  i .  u  )  u  li  r  u  rf  t)  y  <4 
IF  (TIME  (X+  1)  -BB.LT.  0)     G 

9700  BB=3B-24+1C0 

9  6 9 ^  IF  (TIME  (X  +  1)-EB-LE.O)  G 


C  TO  96  J  4 


d9^  IF  (TIME (X  + 1)-EB-LE.O)  GC  10  9b84 
IF  (BB- TIM" 


C 


JIME  (K)  . EQ.O) GC  IC  9  698 
DDI=BB— TIME  (X) 


SB  IT  E  { a  ,  38  1 2 )  B  3 ,  D Dl ,  IT ,  X 
BRITE  (3,3812) B3,DD1, IT,  K 


S812)B3,DD1,IT,  K 
•   F  E  =  »  ,  F 1  0  .  a ,  • 
1  »   K=',I2) 

CALL  PHIDLL  (D 


3*  12  FORMAT  (/, 


DDT=  ',F10.4,«    IT=  ',! 


1  »   K=',I2) 
CALL  PHIDLL  (DD1 , N , L , A , E ,PHI , DEL,D 1 , 02, ND , MD ,1 D) 
CALL  PROD  (PHI,  XKt.,^,  N  ,  1  , XKKH6 , SD , MD , ID) 
h  5 1 1  E  ( 6  ,  4  3  1 2 ) 
4312  FORMAT  (/,  '  XKKM6 


FORMAT  (/,  » 
*  3  3   1 1  £  E  ( K )  • ) 

iRITE 


3  3   1 1 M  E  ( 8 )  * ) 
WRITE  (8,90)  (XKKMb  (J)  ,J=1#S  )  ,BB,  I.IHS  (K) 
WBIT£(4,4812) 

W 8 ITS  (4,30)  (XKKM6  (J)  ,  J=  1  ,  N)  , 35, TIME  (K) 
tKM6(KKK)=XKKM6<1) 

v  v  a  C-.   i  vt  y  v  \  —  v  v  v  \it  C.  i  ';  \ 


XKM6  (KKK)  =XKKM6  (2) 
9  698  XKKK=KKK 


TI  =  I 


X  K  M  6 


K  K  K 


WRITE{4,8969) 
3969  FORMAT (/, '   YKM6 
1   I  ') 

Si  B  IT  E  (  4  ,  9  0 )  ¥  K  M  6  ( K  K  K )  ,  X  K  M  6  ( K  K  K )  ,  X  K  K  F  ,  E  3  ,  T 

KKK=KKK+1 

BB=3 3+6 
9634  CONTINUE 

WRITE (8, 48 12) 

WRITE (3, 90)  (XKXM6 <J)  ,  J= 1 ,N )  ,BB, IIME  (K 


BB 


:k) 


64 


322 


WRITE  (4,  '43  1z) 

KaiTE{4,90)  (XXK«6(J)  ,J=1,N)  ,  BB,II21S  (K) 

IT=  ( II  3E  ( K  + 1 )  -  1 50 0 )  *  ( .  0  1 ) 

IS= (TIME (K +  2) -1500) *  (.01) 


WHITE (4, 9b 9b)  IT 
9696  FORMAT  (5X,(   IT  =  ',14) 

1=24*11 


C 

c 
c 

C   LDI 

C 

c 
c 
c 
c 

C7  134 
C 

c 


3=24*15 

1=TI3E  (K+1) -11*100-1500+1 

S=TIHE  [K  +  2)-IS*10  0-150Q  +  S 

D  T  =  S  -  T 

J  K  =  6 

*iRIT£  (8,7865)  J  :< 

WRITI(4,555) 

WRITE  (8,555) 

DO  3  0  22  1=1,  N 

WRITE  (3 , 90)  (PKKH1  (I , J)  ,  0=  1  , N) 

..RITE  (4,  90)  (EKKR1  {I, J)  ,J=1,N) 
CALL  PHIDEi  (DT,N,L,A,E,PHI,DEL,D1#D2,  MDrJJD,LD) 
CALL  TRANS  {rEl  ,N,  1  ,  EELT,  SB  ,  ME) 
CALL  PROD  (EEL,  DELI,  N,1,N,Q  ,  ND,MD, LD) 
CALL  CCH SI  (9  (1, 1)  ,C#N,N,C,ND,MB) 
,ID.LDLDLELELriDIDLBIDLDLEDLDIDLDIDtELDL3IDIDLDLDI 

WHITE  (h, 544) 

HEITE  (8,54  4) 

DC  7  134  1=1, M 

WSITE(4,90)   (Q  (1,0)  ,J=1,N) 

3RITE(8,90)   (C  d,0)  ,J=1,N) 

***  HAS  STEP  OF  FILIER  REACHED  FOINI  25? 
If  (  (K-25)  .NE.O) GOTO  2222 
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1961 

19  3  2 
C 

C 


***     5100TINE    10     PLACE     ELLIPSE    DATA    IN     FILE 

THE1=.50*ATAN  [2*PKKftM  (1,2)  /  (PKO  1(1,1)  -PKKM1  (2, 

S IG  2  X=  (PKK  a  1  (1,1)  «-PKKM1(2,2))/2.+PKK*!l  (1,2) 
1       /SIN  (2.*TEE1) 

SIG2Y=  (PKKH1  (1,1)  +PKK81  (2,  2)  )/2.-PKKM  1  (1,2) 
1   /SIN  (2. ♦THE 1) 

SX=(  {SIG2X}**. 5) *5 

SY=  (  (SI32Y)  **.  5)  *5 

P 1=3.1 4  159  265/12 

CT=CCS  (THE  1) 

S  I  =  S  I  8  ( 1 3  E  1 ) 

DC  1981  IELLIP=1/25 

XI=IELLIF 

XP  (IELLIP)  =SX*COS  (PI*  XI)  *CT-SY*SI8  (PT*XI)  *ST+XH 

i?  (IELLIP)  =SX*C03  (PI*XI)  *ST+SY*SIN  (*I*XI)  *CT+XH 

S RITE  (19,1982)  XP  (IELLIP)  ,YP  (IELLIP) 

FORMAT  (2F15. 7) 

***  END  OF  ELLIPSE  CALC01ATICN 


))) 


(25) 

(25) 


2222  CONTINUE 
C 

C       ****  R5IIE  IK  11  STUFF 
KF.IIE(8,2224) 

2224  FORMAT  (//,  '  1KH72  (JJ)  ,  IKM48(JJ),  IKM24{JJ) 
1   ,  TKM6(JJ)  ») 

3 RITE  (8,22  23)  (TKM7  2  (J  J)  ,TKM43  (J J)  ,IKM2  4  (J J)  ,TKM6  (J  J) 
1  ,JJ=1,NN) 
WRITE  (8,2225) 

2225  FORMAT  {//,'  TKP72  (MM) f  TKF43(8H),  IKP24  (MM) , 
1   T  K  P  6  ( M  i )  •  ) 

WRITE  (8,2223)  (TKP72  (MM)  ,TKP48  (MM)  ,1KP24  (Mil)  ,TK?6  (MM) 
1  ,  M  M  =  1  ,  N  N ) 
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2  2  <:  3   F  C  E  3  A  T  ( 4  1 1  0 ) 

c 

C      K5ITE(4,7777) 

SHITS  (8,39  3)   (TIHEd)  ,LA"I(I)  ,LONG(I),YH.{I)  ,XH(I)  ,1  =  1 
1   ,NZ) 

WRIT  2  (o,39  3)   (BT1HE  (I)  ,  ELAT  (I)  ,  ELONG(I)  ,YKfl6  (I) 
1  ,XK36  (I)  ,1=1,  NZ) 
39  3  FC8MAT  (1X/5.S  6.2) 
C      r* RITE  (4,77  77) 

RRITE(4,393)   [THE  (I)  ,LAI  (I)  ,L0NG  (I)  r  YH  (I)  ,XH  (I) 
1    ,1=1, NZ) 

KHITE  (4,39  3)      (BTItlE  (I)  ,  ELAT  (I)  ,  3LCJG  (I)  ,YKii6  (I) 
1       rXKfl6(I)  ,1=1, NZ) 
C 

INTIME=INT  (TIHE  (K)  ) 
RLTIME=FLCA1  (INIII1E) 
IF  (HLTIHE. HE. II21E (K)  )  JO^C    399 

CALL    TERR  (ELAI , BLONG ,TI8£, BTItfE ,XB, IB, LAI,  LCNG) 
3  99       CONTINUE 

c     ***  ;;eite  foei  and  ecb2  tc  file  •  ecs  data1 

8EITE(6,32S) 
329   FORMAT  {//,5X, '    JULIAN  HE  EOE 1 

1   E0B2f) 

HR IT E  (6,3  2  8)  (KJULHR  (IEOr)  ,  E0B1  (I  LCD)  ,FCB2  (I5CE) 
1  ,1203=1, NZ) 

WRITE  (16,3  28)  (KJUIHB  (IEOEJ  ,  E0B1  (IECB)  ,E0B2  (IECB) 
1  , IECB=1,NZ) 
323   FCRMAI  (I15,2F15.2) 
C 

DO  3^7  1=1, NN 

ey(1)  =  yk:i6(I)-blai(I) 

EX  (I) =  XKtfo  (I) -BLCNG  (I) 
TI  HE  (I)  =FLCAI  (I) 
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c 

c 
c 
c 

c 

c 
c 


c 
c 
c 
c 

c 


ill    CONTINUE 

WRITE  (4, 39  3)   (BTIME(I)  ,ELAT(I)  ,BLGNG(I)  ,SX(1)  »EX  (I) 

1   ,i=i,kz) 

w RITE  ( 8 , 9 8 58)  (G  1 1  (K)  , K=  1 , 1  0) 
BEIT E  (8 , 98S8)  (TIM E  (X)  ,  K= 1 ,  1  0) 

98  98  FOB  MAT  (!•  10.  4) 
WHITE  (8,410) 
4  10  FORMAT  ( '  1  * ) 

fcRITE(4,7777) 

CALL  PLOTP  (TIM  E  ,  G 1  1 ,  ti  N ,  1 ) 

CALL  PLC  IF  (1IME,G3  1,NN,3) 

WHITE  (8,410) 

WRIT E(4, 7777) 

CALL  PLOTS  (TIME, YH,NN, 1) 

CALL  PLOTP  (TIME, YHF,NN, 3) 

L  L  =  I  IN  +  2 

LONG  (LL-1) =100 

LONG(LL)  =  160 

LAI  (LL-1)  =0 

LAT  (LL) =t0 

WRITE  (4,7777) 

CALL  PLOil  (LONG, LAT, LL,  1) 

CALL  PLC IT  (XH,XH,NN,3) 

CALL  PIOTP  (LONG, LAT, LL,  1) 

CALL  ?LCI?  (XH,IH,NN,3) 


C  *****  WRITE  VAIOFS  INTC  FLCT  FILE  FOR  DI53PLA  ***  ***  **  ** 


:  NUMBER  CF  VALUES 

K E IT E  ( 1 0 ,  2 4 2 8)  N N  ,  K  ( 1  ,  1 )  , E  (1,1) 
WRITE  (  4,2  42  a)   NN,W(1,1),  8(1,1) 

2  42  0  FORMAT  (14, 3X,  ' W  (1 , 1) =» ,F8.3,3X, ' 3  (1  ,  1) 

:  COLUMN  HEADINGS 


=',F3.2) 
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WEITE(10,2429) 


242°    FCRttAl     ('  INDEX  »  ,  3X,  *TiaE  '  ,  41,  'G1  1  *  ,8X,  '  3  j1  •  ,  o ;( ,  '  XH' 
*,^,     •  Y  H  ?  •  ,  5  X  ,  1  L  C  N  3  •  ,  5  X  ,  '  L  A  I  •  ,  7  \  ,  »  X  n  •  ) 
C 
C  VALUES    TO    EE    PIC ITED 

WHITE     (10,2430)  (I,  II HE  (I)  ,G1  1  {I)  ,G31  (I)  ,YB<I)  ,  YHP(I)  , 
*LCNG  (I)  ,LA1(I)  ,XH  (I)  ,    I    =     1,NN) 
C 

2430    FOEHAI     (I4,3F9.2) 
C 

KBITE(8,1016)K 
C  UJtfTE(9,1016)K 

DO    b5    J-1,N 

tfRITE(8,90)   (G(J,I)  ,1=1,-1) 
b5    CONTINUE 

WB ITS  (8, 16) 
1b  FORMAT ('  »,/,3X, ♦CCV.  MAI.  OF  PREDICTED  ESTIMATE'  ) 
CO  1 j  J=1,N 

SPITE (8, 90)  (PKKH1  (I,J)  ,1=1, N) 
13  CONTINUE 
8dd  SIC? 
END 

c  h  h  u  m  h  .i  w  a  s  .i  *  a  h  a  a  w  a » a  w  a  h  a  w  a  a  w  a  w  a  w  a  w  w  a  w  a  u  a  w  .1 .1  *  s  u  m  .•;  .1  *  a  w  a  w  w  a s  a 
c 

SUBROUTINE  PHIDEI  (I,N,a, A, 3, PHI, DEL, D1 , D 2, ND, HI, LD) 
DiaE »S 10 S  A (12,12)  ,8(12,12)  , Pfl I  ( 1 2 ,  1 2)  , DEL  ( 1 2 , 1 2)  , 

*  IE ha  (12,  12)  , 

*  COR  (12, 12)  ,C  (12,12)  ,D1  (12,12)  ,D2  ( 12,  12)  ,TEIL  (12,  12) 

TEST  =  1.E-7 

F=1. 

DO  10  IB  =  1,N 

DO  10  IC  =  1,N 

P  H I  ( I R  ,  I C )  =  0  , 

PHI  (IE f IE)  =  1. 
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C  (IE,IC)   =  h   (IB,IC) 

13II(IH,IC)  =  I/2.00*PHI(IB#IC) 

10  IEBM  (IR,IC)  =  I*EHI(IB,IC) 
5  0    DO  11  IB  *  1,11 

DO  11  IC  =  1,N 
COR<IR,IC)  =  T/E*C(IB,IC) 
PHI<IR,IC)  =  PHI  (IE,IC)  +CC?.  (IB,IC) 
TEI.L(IR,IC)  =  TEIL(IR,IC) +T/(  (F+1.   )*(F  +  2-   )) 
1      *CCE  (IE,IC) 

11  IEP.il  (IB,  IC)  =  TEBDi  (IR,IC)  +1/  (F  +  1.   }  *CCF,  (IB  ,  IC) 
DO  12  IP  =  1,8 

DO  12  IC  =  1  ,  N 
C  (115,1c)  =  0. 
DO  12  K  =  1,N 

12  C(IS,IC)  =  C (IR,IC) +A(IR,K) *COfi  (K,IC) 
F  =  F+1. 

DO  13  1&    =     1,8 
DO  13  IC  =  1  ,  N 

IF (A3S (CGE (IB,IC)  )  .GT.TEST*ABS  (PHI (IR,IC) ) ) 
1     10  IC  5  0 

13  CONTINUE 

CALL  P  R  C  B  { T  E  R  S  ,  B  ,  N  ,  N  ,  K  ,  D  E 1 ,  N  D  ,  D!  D  ,  1 13 ) 

CALL  PROB  (TEILrB,N, H,a# D2,flD,MD, ID) 

DO  14  IB  =  1,N 

DO  14  IC  =  1,  M 
n    D1(IE,IC)  =  DEL(IR,IC) -D2  (IRfIC) 

R  £  T  0  a  8 

5ND 
C  THIS  SUBROUTINE  COMPOTES  IHE  OPTIMUM  GAIN  SAIBIi  AND  1H2 
C  CO VARIANCE 
C 
C 
C 
C 

SUBROUTINE  G  AI8  (PKK,  FKKM  1  ,  Q,R, PHI,  fl,  N  ,  N  ,  G,HI,ND 
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1    ,MD,LD,K) 

DIMENSION  FKK  {  12,  12)  ,Q  (12,  12)  ,H  (12,  12)  ,  J  (12,  12) 
1  ,B  (12,12)  , 

1   [JI<12,12)  ,EI  (12,12)  ,IE8P(12,12)  ,IEMP1<12,12)  , 
1  TEMP2  (12.  12)  , 

S P HI  { 1 2  ,  1 2)  ,  F HI T  ( 1  2  ,  1  2 )  ,  P  '< K  1 1  (12,12) 
C      S  (K)  =  P  (K/K-1)  *HI*  {H*P  (K/K-1)  *ET  +  3) 
C        PHI*P  (K-1/K-1)  *PHI1  +  G 
IF  (K.EQ. 1)  GC  TO  8369 

CALL  T3ANS  (PHI, N , N , PHIL, ND,MD) 
CALL  PROD (PKK,?HIT, N , fc , N , TEMP, NE , VC, I D) 
CALL  PROE (PHI/TEaP,N,K,N/IEap1,NC,aD,LD) 
CALL  ADD  (TEMP  1,  2,N,N  ,  JKK.11  ,ND,MD) 
WRITE  (8, 555) 
b^5    FORMAT  {/'     MATRIX  P8KM1  • ) 
DC  3022  1=1,  N 

WRITE  {8,  90)  (PKKM1  (I,  J)  ,J=1  ,\T) 
3  0  2  2  W BITE  ( 4 , 90 j  (  ? K  K  M  1  ( I ,  J )  ,  J=  1 ,  N ) 
dS'13    CONTINUE 

CALL  TRANS  (H,fl, N, HI, ND,MD) 
C      ERIIE(8,39) 

39  FORMAT  <«      H    •) 
C      CC  22  1=1, fl 
C   2  2  WRITE  (8,90)  (H (I, J)  ,J=1,N) 

9  0  FORMAT  (1P7E1 1.4) 
C      tfEITE{8,36) 

jo  FORMAT  (  ■       HI   » ) 
C      DC  23  1=  1  ,  N 
C  23    KaITE(8,90)  (HT  (I,J) ,J=1,3) 

CALL  PSGD (PKKM1,HT,N,N,M,TEMP,NB,«D,IE) 
CALL  ?  U  0  D  (  H  ,  1  £  .vj  P  ,  H  ,  N  ,  M  ,  T  E  21  P  1  ,  N  0  ,  fi  D  ,  I  D ) 
C      SHITE (8, 38) 

3  8  FORMAT  ('   H  7  H I  '  ) 

C      D  C  2  4  1=  1 ,  M 
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C       24    tfEITE(3,90)  (TEAP1  (I, J)  ,0=1, M) 

CALL    ADD  (TZ/i?  1,?,:1,J#  CE^r 1 , ND, 3D) 

CALL  3ECIP  (*# 0.00  0 COO  1 , TE3P  1  , IEM?2,K£H, HD) 

C      EE.MP2  (1,  1)=1"E;i?3  (2,2) /DEI 

C      IEMP2(2,1) =-lEMP3 (2,1)/CEI 

C      IE3P2(1,2) =-lEHP3 (2,1) /DEI 

C      TE3P2  {2,  1)  =-TEBP3  (1,  1)/BEI 

C      8RIT£(8#31) 

C      DC  2  7  1=1, A 

C   /7  WHITE (8, 90)  (TEaP2(I,J)  ,J=1,X) 
3  1  FORMAT  {'      (HPH  *E)-1») 
I?  (KEE-2)  101,110,101 
1  10  WHITE  {3, 111) 
111  FORMAT  (5HKEF=2) 

101  CALL  PROD  (IEMP/IEMP2,  N,M ,M , G, ND, MD, ID) 
C   NOTE  HEF.E  PKK(I,J)  =  r(K/K)  SHERE  P  (K/K)  = 
C    (I-G  (K)  *H)  *?  (K/K-1) 

CALL  PBCD  <G,H, N  , M, N, IEHE, ND, 8D, IB) 
C      KBITE<3,30) 

C   30  FORMAT  {'     GLi    ') 

C      DC  25  1=1  , N 

C   2d  SBITE  (3,  9  0)  (IE MP  (L  ,J)  ,  0  =  1,  N) 

EC  108  I  =  1,N 

DC  108  J  =  1,N 
Ud  IE -IP  (I,  J)  =  -TEMP {I, J) 
C      WBITE(3,37) 

3  7  FORMAT  ('     HI    • ) 
C      DG  45  1=1, N 
C   4  5  H BIT E  (3 , 9 0 )  (   3 1  { I  ,  J )  ,  J  =  1 ,  N ) 


C 


CALL  ADD  (HI,TEiHP,N#N,T£flE,Ni>,MD) 
WRITE  (8,33) 
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c 
c 
c 

c 


C 

c 

c 


j 3  FORMAT  ('   I  -GH   '  ) 

DC  3  5  I =  1  , 8 
35  WRITE{6,90)  (1E3P(I,J)  ,J=1,8) 

NOTE  HExiE  PKKal(I,J)  =  P(K/K-1)  HfiERE  P(K/K-1) 
CALL  PROD  (IE EP , PKK  H  1 , 8 , N , N , PKK , N C , 3B , I D) 
3  E 1 U  E  N 
2ND 


152 


152 


SUBROUTINE    ADD  (A f E, 8, H,C, 8D,ME) 
DIilE8SrOiN    ft  (8B,i!lB)  ,3  (NB,.MD)  ,C  (2)D,aD) 
DO  2    1=1,8 

DO  2    J    =    1rM 

C  <  I ,  J )     =  0  . 

DO     152    I    =     1,8 
DO    152    J    =     1fH 
C(I,J)     =    A  (I,  J)     ♦    E  (I, J) 
RETURN 

END 


SUBROUTINE  3UB(ArB,8,MfC, »D,MD) 

DI.MEftSIOB  A  (N'D,  IE)  ,  B  (KB,  SB)  rC  (ND,/"C) 

DO  1 52  I  =  1,8 

DO  152  J  =  1,M 

C  (I,  J)  =  A  (I,  J)  -  B  (I, J) 

RETURN 

END 
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151 


1  53 


3 U  BhCU I I N  R    F  RO D  ( A , fi , N , M , L , C , M D / MC , LD) 
DI  *1ENSIC  S    A  (  N  D  ,  M  D)  ,  3  (fi C  ,  L D)  , C  ( ND ,1D) 
DO    1    I    =    1 ,  N 
DC     1     J    =     1,1 

c (i, j)  =  o. 

DO  1d1  I  =  1,8 

DO  151  J  =  1,1 

DC  151  K  =  i,a 

C(I,J)  =C  (1,0)  +  A  (I,K)*B{K,J) 

R  S  T  U  E  N 

END 


SUBROUTINE  TEAMS  ( A,  N,  H,C  ,  ND,flD) 

DI3ENSIC N  A  ( ND , H D)  ,  C  (£D  ,  ND) 

DO  153  i  =  1,N 

DO  153  J  =  1,3 

C(J,I)  -  A  (I,  J) 

BETOSN 

END 


SUBROUTINE  CONST  (Q,ArN,a,C,ND,aD; 
DIMENSION  A  (ND,fiD)  ,C  (KD#HDJ 
IF(Q)   11,10,11 


100 


DO  100  I  =  1,N 

do  100  j  =  i,:i 

C  ( I , J )  =  0.0 
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?  s  i f  j  f  n 

11  IF  (Q-1.C)  13,12,13 

12  DO  120  I  =  1,N 
DO  120  J  =  1,1 

120    C  (I,  J)  =  A  (I,  J) 
B  E  I  'J  F  N 

13  IF  (2+1-0)   15,14,15 

14  DO  140  I  =  1,N 
DO  140  J  =  1,1 

14  0    C  ( I , J )  =  - A  { I , J ) 

EEIOSN 

15  DO  150  I  =  1,U 
DO  150  J  =  1,1 

15  0    C  ( I , J )  =  C  *  A  ( I , J ) 

R  E 1 U  fi  N 
UNO 


SU3BCUTI NE  5 p CI P  {  N , I 9 , A , X  f  K  E  P  ,  1 ) 
DIMENSION   A  (.1,1)  ,X  (1,1) 
DO  1  I  =  1,1 
DO  1  J  =  1,1 

1  X  ( I ,  J )  =  0  . 
DO  2    K    -     1  ,  N 

2  X(K,K)  -  1. 

10  DO   j 4  L  =  1,J 
K?  =  0 

Z  =  0. 

DO  12  K  =  1,N 

IF  (Z.GE.ABS  (&{K,J) ) )  50  IC  12 

11  u    =  AES  (A  (K,  L)  ) 
KP  =  K 

12  CONTINUE 
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1?     (L.GE.KP)  Go  IC  20 

13  DO  1 4  J  =  I ,  S 
L    =    h(Lr  J) 
A(L,J)  =  A  (KF,J) 

14  A (K?,J)   =2 
DO  15  J  =  1 ,  N 
Z  =  X(L,J) 
X(L,J)  =  X(KP,J) 

X (KP,J)  =  Z 

IF  {AE3  (A  (1,1)  ) -LE.EP)  GG  TO  50 
IF{L.3E.U)  GO  TO  34 
L  ?  1  =  L  +  1 
DO  3  6  K  =  LE1,N 
IF  (A  (K,I) •EC.O.)  GG  TC  3o 
BATIO  =  A  (K,L)/A  (L,L) 
DO  3  3  3    =    IP1,N 
A(K,J)  =  A  <K,J)  -EA1I0*A  (L,J) 
DO  3  5  J  =  1,N 

X  (K,  J)  =  X  {K,  J)-BAIIO*X  (L,  J) 
CONTINUE 
CONTINUE 
DO  4  3  I  =  1,M 
11  =  N+1-I 
DO    4  3  J  =  1,N 
5  =  0. 

IF(II-GE.H)  GO  IC  43 
4  1     IIP1  =  11  +  1 

DO  4  2  K  =  II?1,N 

4  2    3  =  S  +  A(I1,K)*X(K,J) 

43    X(I1,J)  =  (X  (11,  J)-S)/A{I1,I1) 
3  E I U  B  N 

5  0    KEP  =  2 

B3T0EH 
E  N  D 


15 
20 
3  0 
31 


32 

33 

J5 
36 

34 
4  0 
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c 
c 
c 
c 


10 
zO 


10 
2  0 


SUBRODTIilf     vlnEAD  (A  ,  N  ,  M  rN  D  ,  DID  ,  IRE  AD) 
DI 3E  3  SIO  S    A  (  N  D , 3 C )  , I f I  A  E  ( 1  0 ) 
DO     10     I    =     1,N 

READ(5,20)  (A  (I, J)  ,J  =  1,M) 
FORMAT  (3E10. 5) 
RETURN 
END 


SUBROUTINE  HKSITE  (A, N , S,ND, 3D, IWBITE) 

DIflENSICIs  A  (ND,MD)  ,1  WHITE  (10) 

DO  10  I  =  1  , N 

«RITE(4,2Q)   (I, J, A  (I, J)  ,  J  =  1,8) 

WHITfi(8,20)   (I,J,A(I,J),  J=  1,8) 

FOR8AT(2(3X,»  (M2,*,»rI2,*)  =  ',1PE10.3)) 

RETURN 

END 


SUBROUTINE  TE~i<  (3  LAI  ,  BLCNG  ,TIKE#  £TIME  ,  Xti,  YH,  LAT 
1  ,LONG) 

DIMENSION  IIME  (30  0) ,BTI8E  (300) ,ELAT  (300) ,BLONG (300) 
1  ,711(300) 

DIMENSION  XH  (3  00)  ,131  (30  0)  ,  EB 2(300)  , IE  JUL  (300) 
1  ,ISJUI(300), 

*  JT.I8E  (300 )  ,  EO  b  1  ( 3  0  J )  ,  EC 3 2  (300) 
PEAL*4  LAT  (300)  , LONG  (300) 
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CALCULATES  rrFCo  121  POSITION  OF  KAL3AN  F.ILTEE 
P3  EDICT  ION  SAND  BEST  I  ft  AC*  V  ALOES  AS  A  FDNC  II CN  OF 
JULIAN  TIM,F  AND  WRITES  THE  DATA  TG  THE  FILE 
•KKER3S  DAIA'. 


C 

C    *** 

C    *** 

c    *  *  * 
c    *** 

c 
c 

DO    26    1=  1, 50 

1 5  T  I  /j  E  =  I  N I  ( 1 1 3  S  ( I )  ) 

CALL    JTJTIHE  (ISTIS E  ,  JUIHE) 

IS JUL  (I)  = JUL  HE 

DO    10    J=1,50 

IBEST=INT (EI IMS (J)  ) 

CALL  JUTIME  (IB2SI,I3J) 

I3JUL  (J)  =IEJ 

IF  (IBJtJL  (J)  . NE-ISJOl  (I)  ) GOTO     10 

EB1  (I)  =  YB  (I)  -BLAI  (J) 

EB2  (I)  =XB  (I)  -BICKG  (J) 
C  2031  (I)  =YH  (I)  -LAI  (I) 

C  EOb2{I)  =XH  (IJ-LCN3  (I) 

EOB1  (1)  =IAT  (I)  -ELAI  (J) 

EOB2  (I)  =  LCNG  (I) -ELONG  (J) 

JTI'IE  (I)  =IEJUL  (J) 
10      CONTINUE 
2b      CONTINUE 

UEIIE  (14,290) 

SBITE  (o,  29  0) 

290  FOR.1AT  (//,  13X,  » JTI21E'  ,5X  ,  »EB1  ',6X,  »f  B2»  ,  5X, 

1        •  SOB  1  » ,5X,  • ECB2' ) 
DO    3  10    N=1,5  0 
IF{JTIMS  (N) .EQ.OJGCTG    210 
WRITE  (6,  200)    (JTIME  (N)  ,EB1  (N)  ,  EB2  (N)  ,20B1  (3) 

1        ,EC3i  (N)  ) 

a  P.  IT  £  ( 1  4  ,  3  0  0 )   ( JT 13 E  (N )  , E 3 1  { N)  , E 32  { N)  , EOB  1  ( N) 
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1   ,E032(N)) 
300     FORMAT  (15X,I9,4F9. 2) 
31)     CONTINUE 

REIT  B N 
FUD 


C 

c 

c 
c 
c 
c 


SUBROUTINE  JUTIME  (ITIBE , JULfcE) 
*********************  LOGNC  E  C  U 1 1 N  E  ********** 

***  CALCULATES  JULIAN  TIME  FPO:-!  YEA?  1900 

***  IYF  =  YEA:\,  iaC  =  5IONTH  (1AECH)  ,  IEA=DAY, 

***  IHR=HCUE  QE  DAY 

I  If  B  =  1 9  6  2 

1,10=3 

IDA=ITIi1E/100 

IHB=ITI3£-IDA* 100 

CALL  NUSCEN  (IYR,IMC#IEA#IHR, JUIHE) 

r<  s  i  o  s  n 


603 


END 

SUBROUTINE  ND3C2  8  (YEAR, BO, DAr HR , JULHR) 

***  CALLED  BY  SUBROUTINE  JOIIHE 

***  CALCULATES  JULIAN  DAY  ANE  JULIAS  HCUF 

IMTEGEE  I N I D  ( 1 2 )  ,  I E  A  R  ,  0  A  ,  H  E 

CAIA  IN  IE/0, 3  1,59,9  0 ,120, 151, 181,212,243,273,3  04 

,334/ 

ID=     (YEAP-1900.) *365.25-0.25 

IADD=0 

IF  (BOE  (YEAB,  4)  .GT.  C)GCTO  b03 

IADD=1 

JULOri  =  INID  (SO)  +  DA  +  LADD 

I H  B  =  24.  *  (10  +  J  U  L  D  A  -  1)  +  HP  +  0.5 

I Y  E  A  B  =  Y  E  A  fi  -  1  9  0  0 
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JULH5=In5 

r  z  i  a  f  :; 

E  N  D 
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